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Die nicht-Riemannsche Geometrie der Versetzungen von KoNpDo sowie von BILBy, 
BULLOUGH und SMITH wird zu einer nicht-linearen Elastizitatstheorie der Versetzungen 
und Eigenspannungen vervollstandigt. Die Grundgleichungen beinhalten formal das 
Verschwinden des zu dem nicht-Riemannschen natiirlichen Zustand des Kérpers mit 
Versetzungen geh6rigen Einstein-Tensors; physikalisch sind sie die Bedingungen fiir 
den Zusammenhalt des Korpers (§ 2). Es gelingt, diese Bedingungen in Spannungs- 
funktionen und Versetzungsdichten allein auszudriicken, und zwar in den Koordinaten 
des deformierten Zustands. Da mit der Einfiihrung der Spannungsfunktionen die 
Gleichgewichtsbedingungen identisch erfiillt sind, reichen die genannten Gleichungen 
zusammen mit den Randbedingungen und dem nicht-linearen Elastizitatsgesetz aus, 
den Spannungszustand des K6rpers aus der Versetzungsdichte zu bestimmen (§ 3). 

Als Beispiele werden die Grundgleichungen fiir Versetzungsanordnungen, die von 
der z-Koordinate nicht abhangen, explizit angegeben (§ 4) und die Spannungsfelder 
der singularen Schrauben- und Stufenversetzung im unendlich ausgedehnten, iso- 
tropen Medium berechnet (§ 5, 6). In der Diskussion (§ 7) wird gezeigt, daB Probleme, 
die sich in der linearen Theorie besonders gut mit Spannungsfunktionen behandeln 
lassen, diese Eigenschaft in der quadratischen Theorie beibehalten. Insbesondere 
ergibt sich eine sehr iibersichtliche Darstellung fiir die mit der Airyschen Spannungs- 
funktion zu behandelnden Probleme. Temperatur-, Konzentrations- und magneto- 
striktive Spannungen lassen sich leicht in die Theorie einordnen. 


Eine allgemeine nicht-lineare Theorie der Eigenspannungen unterscheidet sich 
von der konventionellen nicht-linearen Elastizitatstheorie prinzipiell dadurch, daB 
in ersterer die (nicht-linearen) Kompatibilitatsgleichungen nicht erfiillt sind. 
Infolgedessen wird die in der konventionellen Theorie [/, 2, 3] beniitzte Vorstel- 
lung eines elastischen Verschiebungsfeldes gegenstandslos, und man hat zur 
Bestimmung der Eigenspannungen bei gegebenen Eigenspannungsquellen die 
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Methode der (im allgemeinen dreidimensionalen) Spannungsfunktionen zu be- 
niitzen. Diese Bemerkung gilt nicht fiir den Fall, daB die Eigenspannungsquellen 
im Kérper singular verteilt sind, denn dann existiert auBerhalb der Singularitaten 
ein (wenn auch im allgemeinen mehrdeutiges) Verschiebungsfeld, mit Hilfe dessen 
man nach der konventionellen Theorie arbeiten kann. Anwendungen dieser 
Theorie auf gewisse besonders einfache Versetzungsprobleme sind kiirzlich ge- 
geben worden [4, 4]. 

Sowie man zu etwas weniger einfachen Problemen iibergeht, erweist sich der 
Formalismus der konventionellen Theorie als sehr schwerfallig. Die Erfahrungen 
in der linearen Theorie haben definitiv gezeigt, da8 Versetzungsprobleme sich 
im allgemeinen wesentlich einfacher mit Hilfe der Spannungsfunktionen als mit 
Hilfe des Verschiebungsfeldes behandeln lassen. Man darf daher wohl einige 
Hoffnung an eine Methode der Spannungsfunktionen in der nicht-linearen Elasti- 
zitatstheorie nicht nur stetig verteilter, sondern auch singularer Versetzungen 
kniipfen. 

Eine solche nicht-lineare Theorie der Versetzungen oder — was auf dasselbe 
hinauskommt — der Eigenspannungen kann man in zweierlei Weise entwickeln. 


4, Analog zu der von KRONER [6] entwickelten linearen Theorie der Eigen- 
spannungen, indem man die dortigen Schliisse nun im Hinblick auf groBe Ver- 
formungen wiederholt. Diese Theorie bedient sich einfacher GréBen, namlich 
neben Versetzungen und Spannungsfunktionen im wesentlichen der Deforma- 
tionen, Drehungen (deren Summe als Distorsionen bezeichnet wird) und des 
inzwischen gut verstandenen, fiir den Eigenspannungszustand fundamentalen 
Inkompatibilitatstensors. Die Theorie hat den Stil der Darstellungen der Elasti- 
zitatstheorie in den gebrauchlichen Lehrbiichern, von denen etwa das schéne 
Werk von C. B. BIEZENO und R. GRAMMEL [7] besonders erwahnt werden mége. 
Fur den Kenner dieses Buches sollte die lineare Theorie der Versetzungen leicht 
zu begreifen sein. Die Ausdehnung der Theorie auf groBe Verformungen hat 
vor allem in einer Trennung der sich zunachst fiir die Distorsionen ergebenden 
Gleichungen in solche fiir Deformationen und Drehungen allein zu bestehen. Da 
sich nun Deformation und Drehung nicht mehr linear iiberlagern und die Rota- 
tion kein antisymmetrischer Tensor mehr ist, birgt diese Art von Analyse man- 
cherlei Schwierigkeit. Wir werden daher den zweiten Weg beschreiten, obwohl 
es durchaus méglich ist, da man spater einmal den ersten Weg als den ein- 
facheren ansehen wird. 


2. Der zweite Weg besteht in der Anwendung der nicht-Riemannschen Geo- 
metrie der Versetzungen, wie sie unabhangig von Konpo u. Mitarb. einerseits [8,9] 
und von Bitsy, BuLLouGH und SMITH u. Mitarb. andererseits [10—13] in den 
letzten Jahren entwickelt wurde. Diese Theorie hat den Vorzug, daB sie von vorn- 
herein fiir groBe Verformungen gilt und ist bereits soweit ausgebaut, daB man nur 
sehr wenig hinzufiigen muB, um zu den geometrischen Grundgleichungen fiir die 
Eigenspannungsbestimmung zu kommen. Die nicht-Riemannsche Geometrie der 
Versetzungen ist eine héchst eindrucksvolle und in sich geschlossene Theorie von 
groBer Allgemeinheit. Andererseits wird der mit der Differentialgeometrie weniger 
vertraute, eine Eigenspannungsberechnung beabsichtigende Leser die Beniitzung 
des Formalismus des Ricci-Kalkiils als eine kaum zumutbare Belastung emp- 
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finden. Hierzu ist indessen zu sagen, daB man zwar fiir die Ableitung der geome- 
trischen Grundgleichungen der Eigenspannungsbestimmung die Begriffe der all- 
gemeinen nicht-Riemannschen Geometrie braucht, daB in den Gleichungen selbst 
jedoch nur wie in der linearen Theorie die gebrauchlichen Begriffe Deformation, 
Versetzungsdichte, Inkompatibilitat und Spannungsfunktion vorkommen. In 
den Bezeichnungen halten wir uns hauptsidchlich an das Buch von SCHOUTEN 
([74], insbesondere Kap. III, §1 bis 4). 

In §1 und §2 wird einiges iiber die Idee des nicht-Euklidischen natiirlichen 
Zustands eines Mediums mit Eigenspannungen gesagt. Eine ausfiihrlichere Dar- 
stellung dieser Gesichtspunkte findet man bei Konpo [9]. In § 2 werden ferner 
die genannten geometrischen Grundgleichungen der Eigenspannungsbestimmung 
bei gegebener Versetzungsverteilung abgeleitet. Man vergleiche hierzu insbeson- 
dere die Arbeiten von Konbo und von BixBy u. Mitarb. § 3 hat einiges tiber die 
Lésung der Grundgleichungen im allgemeinen zu sagen, § 4 bringt die Speziali- 
sierung auf ebene Zustande, von denen in §5 und §6 vor allem Anordnungen von 
Schrauben- und Stufenversetzungen untersucht werden. Dabei wird das Problem 
der singularen Schrauben- und Stufenversetzung im unendlich ausgedehnten, 
isotropen Medium im Rahmen der Elastizitatstheorie zweiter Ordnung exakt 
gelost. 

§ 1. Der geometrische Ursprung der Eigenspannungen 

Zur Analyse des Zustands eines mechanisch beanspruchten Mediums hat man 
dieses in zwei Zustanden miteinander zu vergleichen: 

4. Im natiirlichen Zustand. 

2. Im deformierten Zustand. 


Handelt es sich um reine Lastspannungen (auBere Spannungen), so kann 
man den natiirlichen Zustand des gesamten Mediums von dem deformierten aus 
durch Wegnahme der 4uBeren Krafte erreichen. Handelt es sich dagegen um 
Eigenspannungen, so hat man bekanntlich den deformierten K6orper in seine 
Volumenelemente auseinander zu schneiden; erst dann kénnen diese ihren natiir- 
lichen Zustand annehmen. 

Man kann formal das Zustandekommen von Eigenspannungen auf einen 
Zwang zuriickfiihren, der dem Korper das Verlassen des realen (Euklidischen) 
Raumes verbietet. K6nnte man diesen Zwang pl6tzlich beseitigen, so wiirde sich 
der Kérper in einen gewissen nicht-Euklidischen Raum hinein vollkommen ent- 
spannen; man kann m.a.W. auch dem Korper als ganzem einen nattirlichen 
Zustand zuordnen, dieser ist aber nicht mehr Euklidisch. 

Wir kennzeichnen die gegenseitige Lage zweier benachbarter Punkte im defor- 
mierten Zustand durch einen Vektor (Summationskonvention!) 


dy = dx" e,, (1) 


wo die Basis e, ein beliebiges, im allgemeinen krummliniges und nicht-orthogo- 
nales Koordinatensystem (k) des Euklidischen Raumes reprasentiert. Die gegen- 
seitige Lage derselben Materiepunkte im (nicht-Euklidischen) natiirlichen Zu- 
stand sei durch den Vektor dt, = dx*e, (2) 


in einem ebenso willkiirlichen Koordinatensystem (x) des nicht-Euklidischen 
Raumes beschrieben. 
7* 
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Das Quadrat des Abstandes der genannten Punkte ist im deformierten Zustand 


ds? =a, dx dx', api = Cy % 3) 
und im natiirlichen Zustand 
dAst= pdx dx, Cape Bare (4) 


Man kann im letzten Fall die beiden Punkte auch durch ihre gegenseitige Lage im 
deformierten Zustand kennzeichnen und hat dann 


dsp = 85,40" dx', £41 = AR AI bu a- (5) 
Der Zusammenhang zwischen den Zustanden (k) und (x) sei also durch 
dit= Ata, dx a Abd (6) 
gegeben, wobel 
AeA = bP AA Oo (7) 


Wenn Eigenspannungen wirken, sind die Beziehungen (6) nicht mehr wie in der 
nicht-linearen Theorie der Lastspannungen integrabel, d.h. die A haben nicht 
die Form einer Jakobischen Matrix 0 x*/0x*. Ihre physikalische Bedeutung bleibt 
aber offenbar dieselbe: Eine Transformation des natiirlichen Zustands in den 
deformierten Zustand. 

Wir wollen weiterhin alle FeldgréBen als Funktion der Koordinaten x* an- 
sehen. Sinngema8 definieren wir dann den Eulerschen Deformationstensor (é; ;) 
durch die Gleichung 

dst—dst=2e,,dx* dx, (8) 
also 
Sei = Fei — 2Ex1- (9) 
Wir definieren nun die folgenden Ausdriicke 
Paln= 2 Om bey Oommen Creal 
Sintip= @ Vnbkt + Vis —V ebim) - 


Diese Definitionen sollen auch dann gelten, wenn man iiberall g durch a oder ¢ 
ersetzt. Fiir spater brauchen wir ferner die GréBen a”“* und g”*, die wir durch 


(10) 


a”? a, ,= OF, genie OF (11) 
definieren. In den Gln. (10) sei 0,,=0/0x”, wahrend V,, das Symbol fiir ko- 
variante Differentiation im deformierten Zustand beziiglich der a, ,-Metrik ist, 


also Ves Be, tae , 
Vip Shia Cmeenagiten ket eieeres (12) 


Nun ist die Metrik des deformierten Zustands Euklidisch, d.h. der zugehorige 
Riemannsche Kriimmungstensor 


awh is 2[8, 41k — Gl dang Late lnm (13) 


* Wo nicht anders bestimmt, bezieht sich das Heben und Senken der Indizes 
aut die a,1-Metrik. Man beachte, daB diese Operation mit der kovarianten Differen- 
tiation V,, vertauschbar ist, micht aber mit der gewohnlichen Differentiation 6,,. 
Ausgenommen von dieser Vereinbarung sei durchgehend der zu 8pq inverse Tensor oY, 


fir den nach Gl. (11) nicht Vip Ug phd gs sondern i» &jq ghi—g,. gilt. 
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verschwindet identisch. Mit dem an die eckige Klammer angehangten Symbol 
[nm] ist gemeint, daB der antisymmetrische Teil beziiglich der Indizes , m zu 
nehmen ist. 


Die Metrik des natiirlichen Zustands ist hingegen bei Vorliegen von Eigen- 
spannungen grundsatzlich nicht-Euklidisch, d.h. der zugehérige Riemann-Tensor 
Ce 2 [On Smik— Bo" Sonen Snip nan) (14) 
verschwindet in diesem Fall nicht. 
Wegen der bekannten Symmetriebeziehung 


Grnik— = Oni Cont = Gina (15) 
verliert man offenbar nichts, wenn man aus G,,,,;, den zugehdrigen Ricci-Tensor 
Gni= Gumi (16) 


und aus diesem wiederum den zugehérigen Einstein-Tensor 
mi __ goml 1 ml p__ 1 alnt .mjpk 
Gat Gy ae 8 Ga, (17) 


bildet. Die leicht zu verifizierende letzte Beziehung zeigt (vgl. [Zo}s- 454); 
daB 7’ bei kleinen Deformationen mit dem schon frither zur Eigenspannungs- 
bestimmung beniitzten sog. Inkompatibilitatstensor zusammenfallt ([/6], [6], 
S.127). Wie dieser ist 4”’ quellenfrei (V,,”’=0). 

Es empfiehlt sich, in Gl. (17) unter Beniitzung von (412) g,,), anstelle von 
Smlz USW. einzufiihren. Man erhalt so, wenn man gema8 Gl. (9) 


Sale = emt k (18) 


setzt (denn es ist nach Gl. (10) @,,;,=0, da dort nur kovariante Ableitungen der 
4, -Metrik auftreten), 


ni = Sra uel as Emik 1 He Enka Emi] ‘ we 


wobei g?! nach Gl. (9) und (41) der zu (a,;— 26, ;) inverse Tensor ist. 

Setzt man in Gl. (19) 7'7=0, so hat man die bekannten (nicht-linearen) 
Kompatibilitatsgleichungen der Elastizitatstheorie vor sich [2,17]. In einem 
K6érper mit Eigenspannungen ist immer ni 0. 


§ 2. Die Grundgleichungen der Eigenspannungsbestimmung 
Uber den natiirlichen Zustand des Mediums mit Versetzungen wissen wir 
bisher nur, da er nicht-Euklidisch ist und die Metrik (4) hat. Einen nicht- 
Euklidischen Raum pflegt man mit Hilfe des Begriffes der Connexion zu defi- 
nieren [14], und zwar geniigt in unserem Falle die lineare oder affine Connexion 
- ,”(x*) (auch als Affinitat bezeichnet). Mit Hilfe dieser definiert man ein Parallel- 
verschiebungsgesetz fiir Vektoren: Zwei Vektoren ©(P) und €(Q) in zwei um d x” 

auseinander liegenden Punkten P und Q, deren Komponenten sich um 


dC* = —T "Cd xt (20) 


unterscheiden, werden als parallel bezeichnet. 
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Wir gehen jetzt wieder zu dem Koordinatensystem (k) tiber und definieren 


ze . . 
die kovariante Differentiation (Symbol V,,,) beziiglich der Connexion des gleich 
an dem Metriktensor, wobei wir ausnahmsweise jetzt den oberen Index auf die 
&,.-Metrik beziehen, d.h. 

Date Spr = Unk: (21) 


Dann soll gelten 


Vie Sip Operem—Lont Poe ooh ae 

: b b 

Vi Sink= Cio k aR he Cone (22) 
— Vj, 81m = — 81m + Lin? Spm + Dim? &1p- 


Die kovarianten Ableitungen in (22) verschwinden, da sie Ableitungen des eigenen 
Metriktensors darstellen. Bildet man die halbe Summe der Gln. (22), so erhalt 
man leicht, wenn man noch J;,," in symmetrischen Teil toe und antisymmetri- 
schen Teil Ij," aufgespaltet ([74], S. 134ff., [9], S. 462, [18], S. 66) *, 


Liu Smt Past (23) 
mit 


Amtr= Linnea t+ [temi— Liem: (24) 


Den in den Indizes m,/ antisymmetrischen Teil Jf,,;, einer Connexion I}, ;, be- 
zeichnet man nach CarTAN als Torsion [19]. Wir werden sehen, daB sowohl 
Iiniyx als auch h,,;, eine einfache physikalische Bedeutung besitzen. Im Moment 
vermerken wir aber nur, da8 beide GrdBen die Transformationseigenschaften 
eines Tensors dritter Stufe haben, im Gegensatz zu dem symmetrischen Teil 
der Connexion J(,,),, der einer komplizierteren Transformationsformel gehorcht, 
die wir hier aber nicht bendtigen. Man priift weiter leicht nach, daB h,,,;, in 
den letzten beiden Indizes antisymmetrisch ist. 

Der zu der Connexion J}, ,, gehérige Riemannsche Kriimmungstensor schreibt 
sich 


Ripik= 2 Cpl pn eae ele (25) 


und hat die gleichen Symmetrieeigenschaften wie G,,,,;, mit der Ausnahme ([/4], 
S. 138 ff.) 
Ramin He Rika (26) 


Dies bedeutet, daB man wiederum nichts verliert, wenn man zu dem Einstein- 
Tensor tibergeht, daB dieser aber nun im Gegensatz zu nit nicht mehr sym- 
metrisch ist. 

Wir schreiben zunachst den Riemann-Tensor an, indem wir in (25) J7,,, nach 
Gl. (23) und (48) durch ¢,,,, und h,,;, ersetzen. Man erhalt so 


Ry miz= 2[V,(— 2 fry hy 1x) peri 2Enng hake) (=2&mip + lmtp) lbw: (27) 


Bitsy und Smits [12] haben das Verschwinden des Kriimmungstensors R,, nip 
als ein physikalisches Erfordernis bezeichnet, welches notig ist, damit man an 
jedem Punkt eines deformierten Kristalls noch eine Kristallstruktur eindeutig 


* Gl. (23) stellt bekanntlich die allgemeinste metrische Affinitat dar. 
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feststellen kann. Nach der Analyse von KRONER [6] liegt es nahe, das Ver- 
schwinden von R,,,,;, als Zwangsbedingung fiir den Zusammenhang des Korpers 
zu verstehen. Beide Feststellungen meinen denselben Sachverhalt und sind 
makyoskopischer Natur. (In der mikroskopischen Theorie hat man Differenzen- 
gleichungen anstelle von Differentialgleichungen und etwas abgednderte physi- 
kalische Bedeutung, [6], S. 94ff.) Konno hat tiber die genannten Autoren hinaus 
den Fall des nicht-verschwindenden Kriimmungstensors R,,,,;, besprochen ([9], 
S. 465 ff.). Wahrscheinlich hat man damit Kérper zu beschreiben, die makro- 
skopische Risse zeigen. Es wiirde uns hier indessen zu weit fiihren, diese Uber- 
legungen weiter zu verfolgen. 

Wir werden daher das Verschwinden des Riemann-Tensors (27) oder, aqui- 
valent, des zugehGérigen Einstein-Tensors als die geometrischen Grundgleichungen 
unserer nicht-linearen Theorie der Eigenspannungen ansehen. Indem wir sie in 
der Einsteinschen Form anschreiben, erhalten wir fiir den symmetrischen Teil 


$ {el 2 10, (—2 Emin t+ min) —8° 1 ( 2Enngt Muna) (—2Emip + Mmip) 1}a;y=0 (28) 


mit (77) =symmetrischer Teil beziiglich 7,7. Den antisymmetrischen Teil multi- 
plizieren wir mit ¢,;; und erhalten 


t Side Feasae ens V,, Or ae Ns 2 Enkg =i hy nag) on” = Or (29) 


Hier haben wir, vielleicht etwas inkonsequent, aber zweckmaBig, gleich den von 
Konpbo [8] sowie von BiLBy,- BULLOUGH und SMITH [J/] gefundenen wichtigen 
Zusammenhang zwischen Torsion Jj,,;, und (lokaler) Versetzungsdichte «”, be- 
niitzt, in Gleichungsform 


anml 
Dip Snyp Oees an, =— 8” Lint) (30) 
geschrieben, woraus 


ml 
himtk=— % lenip ab, = Eaeh ap ae, — E14 p Com) » ag It (31) 


folgt. Hier ist der (urspriinglich von NYE [20] definierte) Versetzungstensor «”, 
so gewadhlt, daB der erste Index (m) den Linienverlauf, der zweite den Burgers- 
Vektor charakterisiert. 

Nach Gl. (31) kann man h,,;, als eine spezielle Beschreibung der Versetzungs- 
dichte ansehen. Es gibt indessen noch eine direktere Interpretation: Wegen der 
Antisymmetrie von h,,;, in k,/ liegt es nahe, den zugehérigen Tensor zweiter 
Stufe zu bilden, den wir —2K”,, nennen wollen, also (siehe [12], Gln. (19), (20)) 


Ninth = — Enter Kim K", = — 4 hyip: (32) 
Die Beziehungen 


sind dann leicht verifiziert. Der Tensor K”,, ist als Nyescher Kriimmungstensor 
bekannt. Er gibt die Anderung der Kristallorientierung d#” beim Fortschreiten 
um eine Strecke dx” an (und zwar den ,,Korngrenzenanteil"‘, [6], S. 35 ff.), also 


do” = K" dx”. (34) 


Dies folgt aus Gl. (20) (man ersetze J’ durch h). Zu dem Nyeschen Kriimmungs- 
tensor vergleiche insbesondere auch EsHELBY [21]. 
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§ 3. Die Eigenspannungsbestimmung aus den Grundgleichungen 
Der Tensor «*! in den Gln. (28) und (29) gibt die sog. lokale Versetzungsdichte 
an [22, 10, 11]. Sie entsteht aus der natiirlichen Versetzungsdichte des nattirlichen 
Zustands durch die Abbildung auf den deformierten Zustand 


ot? — AP AL gt, (35) 


In der Kristallphysik wird oft eine andere Kennzeichnung der Versetzungen 
gewahlt. Man vergleicht nach FRANK ((22], [23], S. 387) einen geschlossenen 
Umlauf im deformierten Kristall mit einem korrespondierenden Umlauf in einem 
sog. Idealkristall und definiert den SchlieBungsfehler im letzteren Umlauf als 
wahven Burgers-Vektor der Versetzung. Wéahrend eine kontinuierliche Verset- 
zungsdichte wohl im allgemeinen als lokale Dichte vorliegt, ist der Burgers-Vektor 
von singuldren Versetzungen meist als wahrer Burgers-Vektor gegeben. In diesem 
Fall gibt es gewisse Komplikationen, die so gelést werden kénnen, da man mit 
Hilfe der zuerst berechneten linearen Naherung (fiir die der Unterschied wahr — 
lokal noch nicht besteht) die wahre Versetzungsdichte naherungsweise in die 
lokale umrechnet, bevor man den zweiten Naherungsschritt macht usw. 

Bekanntlich kénnen Versetzungslinien im Innern eines K6rpers nicht endigen, 
ohne da8 dort der Zusammenhang zerst6rt wird. Diese Bedingung schreibt sich 
nach NYE [20] im Idealkristall, d.h. fiir die wahre Versetzungsdichte o* 


V0? = 0. (36) 
Fiir eine singulare Versetzung, die mathematisch durch 
ees a (37) 


zu beschreiben ist, (= Tangenteneinheitsvektor der Linie, 6 = (b°) = Burgers- 
Vektor) lautet Gl. (36) wegen V,t*=0- auch 


dbldt =0; (38) 


in Worten: Der Burgers-Vektor einer Versetzungslinie ist langs dieser konstant. 
Ersichtlich ist (29) die Darstellung des Gesetzes (36) im deformierten Zustand 
(fiir die lokale Versetzungsdichte). Fiir eine einzelne Versetzung 7”b* erhalt man 
dann leicht die von der Deformation bewirkte Inkonstanz des Burgers-Vektors 
im deformierten Zustand (des lokalen Burgers-Vektors) zu 


db*|dt = — g°9 (2&4 ng — $ Engr’ On) t” Dy. (39) 


Die weiteren Rechnungen dieser Arbeit fiihren wir einfachheitshalber in 
cartesischen Koordinaten durch. Dann verschwindet der Unterschied zwischen 
Ko- und Kontravarianz, es wird V,,=0,, @jqj,=0, 4;=6,; usw. Alle Tensor- 
indizes schreiben wir von nun an unten. Ferner verzichten wir vorerst auf eine 
Besprechung des Randwertproblems, d.h. wir nehmen ein unendlich ausgedehntes 
Medium an. Die Beriicksichtigung von Randbedingungen stellt jedenfalls keine 
grundsatzliche Schwierigkeit dar (vgl. § 7). 

Wir besprechen jetzt die Lésung der Grundgleichungen (28) unter der An- 


nahme ,,lokale Versetzungsdichte gegeben‘‘. Von dieser verlangen wir, daB sie 
die Gln. (29) befriedigt. 
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Vernachlassigt man in Gl. (28) die nicht-linearen Glieder, und setzt man ¢,,); 
und ,,;, gemaB Gl. (10) und (31) ein, so bleibt der von KR6NER (([24], [6], S. 34) 
zuerst angegebene Zusammenhang zwischen den Deformationen und der Ver- 
setzungsdichte, oder, wenn man will, zwischen dem (linearen) Inkompatibilitats- 
tensor 7; und der Versetzungsdichte 


pe (eae Silk On 0; Emk)(éj) eos (€ jnm 0 WG mG) = ad nt (40) 
ubrig. Die Gln. (28) sind also die nicht-lineare Verallgemeinerung dieser gut 
verstandenen Gleichungen. 

Zur storungstheoretischen Behandlung der Grundgleichungen (28) bei ge- 
gebener lokaler Versetzungsdichte «,, befriedigen wir zunachst die Gleichgewichts- 
bedingungen (bei Abwesenheit von Massenkraften, o;; = Spannungstensor) 


6;6;;=0 (41) 
identisch durch den Spannungsfunktionsansatz ([26], [27], [6], S. 54 ff.) 
O55 = — Ejnm Sith On 1 Ak (42) 
mit der Nebenbedingung 
85 (45 — ce qaegy tar 6:;)/2G =0 (43) 


(G = Schubmodul, vy = Poisson-Zahl, der K6rper ist elastisch homogen und isotrop 
angenommen). Die zugehérigen Deformationen folgen mit Hilfe des zustandigen 
nicht-linearen Elastizitatsgesetzes. Man erhalt so die ¢,,,, ausgedriickt in den 
Spannungsfunktionen. Die sich nach Einsetzen in (28) ergebenden Gleichungen 
lauten 


AA y;;= 26 |n'; Pls Usa ee (“het Pent Qne) 6, (44) 
Hier entspricht 
AA; — 26(n, + mh 6;;) (45) 
dem linearen Anteil von (28) (gleich Null gesetzt ist (45) mit Gl. (40) gleich- 
wertig [25]). P,,; ist definiert durch 
P= (Ejnm Sith On Op Ema)(64)> (46) 


wobei die «f, denjenigen Teil der Deformationen ¢,, darstellen, der aus den 
Spannungen o;; durch Anwendung des nicht-linearen Anteils des Elastizitats- 
gesetzes folgt (§ 4). Die Q;,; sind definiert durch 


O;, = 4 (Ejnm 1a lL neg i mipa;  Lmik= — 2€min + Mmin- (47) 


Die Rechnung ist nun im Prinzip sehr einfach. Man setze zuerst P,;= Q;;=0 
und berechne die lineare Naherung 7%; aus der Gleichung 


5 Mh 5:4) « (48) 


AAYi= 26 (iby + 5 
Danach laBt sich das zugehérige P?-+ Q?; ausrechnen. Der Ausdruck 3+ 44; mit 


AA; =2G|Ph+ 0+ = (Phe + Os) 6] (49) 
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ist eine bessere Naherung, noch besser ist 7?;+-y;+3; mit 


AAy},= 26 [P+ Oy + <2 (Phat Oke) 6:4 (50) 


iy 
usw. Die lineare Naherung ist oft relativ gut zu erhalten. Wenn wir mit 7?;+ 43; 
abbrechen, haben wir das Problem im Rahmen der Elastizitatstheorie zweiter 
Ordnung gelést. Hiermit mu8 man sich wohl im allgemeinen zufrieden geben, 
da schon die fiir eine weitere Rechnung bendétigten Elastizitatskonstanten nicht 
zur Verfiigung stehen [4]. 

Wir bemerken noch, daB es im Rahmen der quadratischen Naherung geniigt, 
g?? in Gl. (28) durch a,,=6,, zu ersetzen, da sich beide GréBen nach Gl. (11) 
und (9) nur um eine GréBe der Ordnung ¢,; unterscheiden. Ferner braucht in 
diesem Fall fiir die ¢,,,, in den nicht-linearen Termen von Gl. (28) bzw. (44) nur 
der mit dem linearen Elastizitatsgesetz erhaltene Anteil eingesetzt zu werden. 
Beide Vereinfachungen sind sehr wesentlich. 

Die lineare Naherung des Spannungsfunktionstensors hat im Fall der sin- 
gularen Versetzung die besonders einfache Form ([27], [6], S. 72) 

Li= ee (Ej a9 rdLi) Qj = Wes — Esra. §;;)/26. (51) 
Hier ist y der Abstand zwischen dem Linienelement dL; der Versetzung und dem 
Aufpunkt. Die Integration erstreckt sich langs der Versetzungslinie. 


§ 4. Ebene Eigenspannungszustande in der Elastizitatstheorie 
zweiter Ordnung 
Zu der expliziten Darstellung der Grundgleichungen fiir die Eigenspannungs- 
bestimmung im Falle ebener Zustande benétigen wir das Elastizitatsgesetz ex- 
plizit, und zwar in der Form ¢,;= €;(0;;). Man gewinnt dieses Gesetz aus der 
zur Murnaghanschen Darstellung ([1], S.64) der Energiedichte e analogen Form 


1 2 1 L+2M 3 ~ 
— 1G; One o 2M o,0 ING ae 
Guay) ae ek 3 I rant UI (52) 
in der ; 
Oy =O, Ae Ongri- One, (53 ) 
Oren, Op. Gn O PRs 
on = y yy yz we 2 x : (53’") 
xy Oyy Oy x Oz% Ox: Ox x 
F Ox x Oxy Ox sz 
On = | Oxy Oyy Oye ’ (S354) 
Ong Oy» Oz 


die skalaren Invarianten des Spannungstensors sind. Die Deformationen schreiben 
sich dann 


a = ¢,,=—~ 1 — — 
€/00z Ex x 2G(1 +9) ee V(Oy, ate Gy.) tg (54') 
air Loi a 2M (oy ‘De O10, x) ae EV Opes Coz = Gee) , 
1 1 
oe de/00,,= 8 = 2G Ov Aa Oy eal On One G08.) (54’’) 


usw. mit zyklischer Vertauschung der Indizes. Zu dem Zusammenhang der nicht- 
linearen Konstanten L, M, N mit den Murnaghanschen J, m, 1 s. Anhang. 
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Wir schreiben zunachst die linearisierten Grundgleichungen (40) mit 0/6z =0 
aleealso 


a Oy b= re dy yz» ae Ge ee) 0, Ay x» apt i (0, yg OD, Gray 2s (55’) 
— 0,(— 8, a sinh Ep = 0G. — CLD) ORS aaa SP 65”) 
Ne oy e+ oy 2) ome dy (op Oi 2 ‘2 Oy Hy y/ 2, 


0 0 Oy == 
a Oe i Oey Eyy + 20,4 Exy = 0 Ke op Key» (55.4) 


wobei 0,,= 6?/0x* usw. Danach machen wir den aus den Gln. (42) folgenden 


,ebenen“ Spannungsfunktionsansatz (F,,=0,,F, ®,=0,® usw.) 


XX 


y 


esi? —-. Pe . y, 
Cipla Oyy =F, Coy Le; (56’) 
6,,=®,, Oy, =— ®,; O,, = O72, (56) 


mit P= —y,, als Airyscher Spannungsfunktion und O= 0, 7, ,—0,4,, als Span- 
nungsfunktion der Torsion (oder Schubspannungsfunktion, [7], [6], S. 53. Die 
Gln. (57’) und (59’) werden zeigen, da8 sich o,, in den uns interessierenden Fallen 
trivial ergibt.) Mit Hilfe des Hookeschen Gesetzes bringen wir die Spannungen 
bzw. Spannungsfunktionen in die Gln. (55) ein. Die trivialen Integrationen dieser 
Gleichungen fiihren wir gleich durch, wobei wir eine unbestimmt bleibende, in 
den Ortskoordinaten lineare Deformation vorerst Null setzen. Es ergibt sich so 
(A p. Cael oy) 


of, =v AF®+ 2G(1 +») &, (57a) 
AND EGA Ce (57) 
AAF® = — i ihe Oye) 


Die Funktionen é, 7, € bestimmen sich héchst einfach aus 


0, F = by, 05 —— Oy, 
2056 SO ne Oh ere Cas DO Cty Oasys (58) 


(ar er Oey =| Y (0, Oy g — Oy Os) - 


Die nun folgende Integration der Gln. (57) liefert die Spannungsfunktionen 
der linearen Naherung ®°, F°® sowie of,. Danach kénnen wir zu dem zweiten 
Schritt des Naherungsverfahrens iibergehen. Die zugehérigen Gln. (49) schreiben 
sich mit 0/0z=0 

a7, =v AF) + 2G(4 +9) (— & + Q’), (59') 

A@' = 2G(d, 62, — 4, 6€, + Q"), (59”’) 


AAP! =— 25 [0,854 2¢,68y— 204588) +9A chet Ort” (Qe + Qyy)], 59") 


wobei ; 
ip = (Sjnm Sire Lang Tniadeizy» — Ltn = — 28min + Monies (60’) 


en Q’ a a OP. Oy y Q' “ahha Qyy Oxy g a Cry: (60) 
6.0" =n 40" = Oe 
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Die «2, folgen definitionsgemaB aus dem nicht-linearen Anteil des Elastizitats- 
gesetzes, angewandt auf oj;. | 

Die Integration des Gleichungssystems (59) liefert Cae OL und i, womit 
gemaB Gln. (56) alle Spannungen 6; ,= 09; + 04; in der quadratischen Naherung 
bestimmt sind. Dies gilt allerdings nur im unendlichen Medium; im endlichen 
Medium hat man noch fiir die Erfiillung der Randbedingungen zu sorgen, woritiber 
in § 7 gesprochen wird. 

Die fiir die Berechnung der Q;, benétigten I;,,, ergeben sich mit Gl. (10) 
und (31) leicht zu 


Lin=0 
0 ae 0 + Sara @ 
i Bis eae ao) Cory Leu = Oy Exe: Lg — Oy Sy ys Loy = Oy Ey y 
& ae 0 a 
pe ae Oy erga Dy = Cpe yiypate Ona, (61 ) 
0 ss 0 oy a 
Migs ce = PR oon a Ch Pa ale hee Dg He, Sey tr Sg byy Key» 
Dean yar Deg a Mp4) (61 b) 
0 pa: 0 
Oe ee Orga Le a Ley Sees (61) 
0 eRe. 0 
Tye — 2.0, by,— Cyy, Ly, = — 20, bys + Oyy, 
= Son Le (61d) 
Dy Cg yy ey 


Hier wurde beriicksichtigt, daB die «2; den linearen Gln. (55) geniigen, ferner in 
den L},;, mit m=-1=-k die aus Gl. (29) in erster Naherung folgenden Beziehungen 
CNP Pe Oy Cy = 0, Oy, Oy + Oy Oy y = 0. 

Die Gln. (61) zerfallen in mehrere Gruppen, die wir jetzt kurz diskutieren. 

a, b) Zu diesen Gruppen gehéren einmal Stufenversetzungen, die in Ebenen 
senkrecht zur z-Richtung verlaufen und ihren Burgers-Vektor in z-Richtung 
haben («,,, %,,). Ist «,,—a,,—0, so verschwinden stets die J},,, der Gruppe b, 
diejenigen der Gruppe a aber dann nicht, wenn Stufenversetzungen in z-Richtung 


vorhanden sind («,,,, %,,). 

c, d) Zu diesen Gruppen gehéren einmal Versetzungsanordnungen mit Linien- 
verlauf und Burgers-Vektor senkrecht zur z-Richtung («,,, Dain Puen hyp ean te 
schwinden diese Versetzungen, so verschwinden auch alle [},,, der Gruppe d, 
diejenigen der Gruppe c aber dann nicht, wenn Schraubenversetzungen in z- 


Richtung vorhanden sind («,,). Damit sind alle neun Versetzungskomponenten 
%,, elngeordnet. 


Von den angefiihrten Problemen sind in der Kristallphysik die weitaus wich- 
tigsten die zu Versetzungen in z-Richtung gehérigen. Wir werden daher weiterhin 
nur noch die beiden Probleme 


,Behiebige Verteilung von Schraubenversetzungen oe uw 
,,Beliebige Verteilung von Stufenversetzungen «,, Coen) ese 2m) he 


behandeln. Dann verschwinden jeweils die meisten der D},in, und auch die Be- 
rechnung der «@, vereinfacht sich wesentlich. 
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§ 5. Schraubenversetzungen 


Im unendlich ausgedehnten Medium sei die lokale Versetzungsdichte «, , (x, ¥) 
gegeben. Als Teillésung der linearen Gln. (57) hat man wegen &=7 =0 sofort 
o:;—0, F°=0. Deshalb ist von den linearen Deformationen ¢, nur ¢°,, ¢°, von 
Null verschieden. Weiter verschwinden alle Q,;,; auBer 


OF = ! by BRE ad I pep i Dine al [(2, ei) (6, ey) fe (0, e%2) (2, €5 2) | ‘ (62) 
Mit e%,=0%,/2G, &,=0},/2G und (56”) erhalten wir 
Oa De Py — Os (63) 


Die Berechnung der ef, nach Gl. (54) ergibt 
ef, = 2M (B29? + @?) _ NG? 
ef, = 2M(@9? + G9?) — N G9? 


&7, = 2M(DY? + @5”) (64) 
ef, =— ND 
ef, = ee, =0. 


Hiernach erhalt man fiir das Gleichungssystem (59) im Falle von Schrauben- 
versetzungen «,,(%, Y) 


0, = AP — 4MG(t +9) (®P + OF"), (65') 
Beans? 65") 
AAP =— 25 [2M(1 +9) A(@Y + Of) — 
= — PE [ame +») AOe? + 05%) + (2N — Fa) (2.3 — B25). 


1 
Als Beispiel fiir die Lésung dieses Gleichungssystems behandeln wir die sin- 
gulare Schraubenversetzung in der z-Achse (“:2= bd(0), ?=2*+ 97, yp =arctg s ; 
Die Spannungsfunktion ®° ergibt sich in diesem Fall zu ([6], S. 75) 


eer, 
Po — Sq me: (66) 
Damit wird aus Gl. (65’”) 
nto Gp? 3 1 
AAF'=4Kiet,  K=— ay [SMU +9) — 24 Ge] (67) 


und hieraus folgt eindeutig (im unendlichen Medium) 
A ae OG (Ore (68) 


(Ein additives Glied const.y hat nicht die erforderliche Rotationssymmetrie.) 


Mit Gl. (68) haben wir bereits 
Oxet Oy = K/0?, (69’) 


ot, = [vK —4MG(1+)]/e. (69’’) 
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Um die iibrigen Spannungskomponenten zu bekommen, integrieren wir Gl. (68), 

wobei wegen der Symmetrie des Problems Fi\=F1(9) sein sollte. Damit wird 

aus Gl. (68) ES oe at F} = K|o? (70) 
Q 


mit der Lésung Fl= K(Ino)?/2. (71) 
(Beachte: Ing ist wegen Aln 0 =6(0) keine Lésung der homogenen Gleichung 
von (68), und ein additiver Summand const. p in (68) bewirkt eine Verletzung 
der integralen Momenten-Gleichgewichtsbedingungen fiir die Zylindermantel- 
flichen.) Die Spannungen ergeben sich gemaB Gl. (56’) zu 


1s =o1,=K(% a In 9 


PM o4 4 
2 2___ 472 \ 
Fi, =0yy=K (+ —{* Ing] (72) 
Q @ 
Ey = Gzy= K—¢ (4 =2In 9), 
oder nach Transformation in Zylinderkoordinaten und gleich fiir o;;=o%;+ 01; 
angeschrieben 4 Ob 
Og 4x ln oes Cy aa Po 
Ooo = K(1 — In)/e? Ona (73) 
0,,=|[vK —4MG(i+7)]/¢? Ong = O- 


Damit ist der Zustand des unendlich ausgedehnten Mediums mit einer Schrauben- 
versetzung in der z-Achse im Rahmen der Elastizitatstheorie zweiter Ordnung 
exakt und vollstandig bestimmt. 

Schneidet man langs zweier Ebenen z= -+-a (a =beliebige Konstante) auf, 
so hat man dort Krafte anzubringen, wenn keine zusatzlichen Deformationen 
auftreten sollen. Uber eine der Flaichen, z.B. z=-+a, gemittelt verschwindet 
weder die Gesamtkraft (wegen o,,) noch das Gesamtdrehmoment (wegen o,.). 
Solche Erscheinungen sind bei ebenen elastischen Zustanden wohlbekannt. Physi- 
kalisch kommen nur endliche Kérper vor, dann werden diese Unstimmigkeiten 
durch Erfiillung der Randbedingungen beseitigt (vgl. ESHELBY [29]). 


§ 6. Stufenversetzungen 


Im unendlich ausgedehnten Medium sei die lokale Versetzungsdichte «, , (x, y), 
a,4(%, y) gegeben. Als Teillésung der linearen Gln. (57) hat man wegen €=0 
sofort «9,=¢69,=0. Daher verschwinden in (61) alle [},,, auBer denen der 


Gruppe a. Das sind gerade diejenigen, bei denen der Index z nicht vorkommt, 
deshalb bleibt gemaB Gl. (60) von den jj; nur 


Org teliyplege lie bby Seer! Be) BSS) ead bee (74) 
iibrig. In den Spannungsfunktionen geschrieben lautet dies 

1 
Or 32 {4G(4 —29) (a,, 4) —a,, AFP) +2 (4 —9) (4 —2y) | (4F?)2+ (AF})?|— 


75 
—A(Fr. Fy — Ey) + AR, + By AR’, — 2F2, APS}. a 
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Die Berechnung der ef, nach Gl. (54) ergibt 
e¢2—=[L(1+)?4+2M] (AF°)?—2M(F2, F°,—Fo?) +[—2M(1+»)+N0] Fo, AF 


eyy=[L(1+9)*+2M] (AF°)*—2M (Fp, Fy —Fe?) +[—2M (140) +N] Fo, AF 


ef, =[L(1 +9)? +2M P*] (AF*)?— (2M —N) (F?, Fyy—F2?) (76) 
ef, =[(—2M(1+»)+Nv] Fo, AF? 
62, = bp ie O 


und der in Gl. (59’’’) benétigte Ausdruck 


R= Oyy 4+ Onn Evy — 20,62, + vA ec, (77) 
berechnet sich zu 


P=[L(1 +»)8 + 2M(1 + »°)] A(AF%)2— 1a 
— [2M +») — No] [A( Fo, FS,— Fo?) + F°, AF®, + F°, AF®,— 2F°, AP?,). 


Hiernach erhalten wir fiir das Gleichungssystem (59) 


0;,= v AF! — 2G(1 +») {(L(4 +-¥)* + 2M] (AF)? — 


79’ 
— (2M —N) (Ey Eo) ae 
ADL Oy (79’") 
AAF = — 2° {C, (4, AF? — yy AF2) + Cy [(4F2)* + (AFS)*] + se 
wey 

BoC AEP Tyas 8), (FR, AR EARS oF ARS) CA (Amo) 

mit 
1—2 {— 1—2 

Cys amit = ( ae Ge = OMA fy) ON 418 G, ih 


Cy=—2M(1+7)4+Nvr41/8G2, C,=L(1+)?+2M(1 +7). 


Als Beispiel fiir die Lésung dieses Gleichungssystems behandeln wir die singulare 
Stufenversetzung in der z-Achse mit Burgers-Vektor in x-Richtung, also «,,= 
bd(e@). Die Airysche Spannungsfunktion der linearen Naherung ist dann nach 
KoEHLER ([30], [6], S. 74) 
Gb 
0 = 

FPo=Aylnga, A= Cea (81) 
Damit wird aus Gl. (79’’’) (das Glied mit «,, in (79’) gibt zu F? keinen Beitrag, 
s. Anhang) 


AAFA=8(B +2C “a We 
B= 22 sea + 22 +9)% + 2M(1 +9) (1— 29 + 29) +No}, (82) 
G A2 
ee lre 


Die maBgebende Lésung ist 
Fi=(Bing +C a 


“Jing, (83) 


und die Airysche Spannungsfunktion einer Stufenversetzung in der z-Achse mit 
Burgers-Vektor (Betrag }) in x-Richtung wird in der quadratischen Naherung 


0? 
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Hijerzu berechnet man nach Gl. (56’), (57’) und (79’) leicht die Komponenten 
des Spannungstensors Zu 


4 
A see {eproene oe ! 
4% 02 0 0? Q 
—(B+2C) eae VERBS eon o} 
Q 
PRIME a ea: [Gat A(R SCE * (28 SO) ee 
Oyy 0? 7 0? o? | @ 
Best 6Cs19* (8-20) " 91mg} 
Q 
oO Vee he ay {B= 3°) PINE IOI = (85) 
ay 0? 02 o4 0? 
(BA 20) Ht (B+2C) ¥? -2In gh 
Q 


6,,=2Aryle?4 [D | <a (E 4C vin o)|/¢° 
D=2{By—2G A2(1 +) LU +7)? + 2M v*)}} 
E =2G6 A2(1+»){2[(L(1 +»)? +2Mr7] -2M+4+ N} 
Cn Onan: 
In Zylinderkoordinaten ist 
F=(Aosing + Blne + Ccos2¢q) Ine. (86) 


Hierzu berechnet man mit 


Se= aloo + he Opp=E,, %9=—~F, +k (87) 


ee@ ro Pe ep ~~ e@? ep 0 0? vo 
leicht die Komponenten des Spannungstensors: 
0. = A(sin g)/o + [2B lne + C(cos 2g) (1— 41nQ) |/0? 
Op, = A(sin g)/o + [2B(1 — Ing) — C cos 29] /¢? 
Sap = — A(cos gle + 2C “222 (1 — Ing) (88) 
0,,—= 2A (sin g)/e + [D + (cos 29) (E — 4C vIn @)|/¢? 


C5, = On, — OF 


§ 7. Diskussion 


In unserer Diskussion der nicht-linearen Elastizitatstheorie wollen wir zu- 
nachst darauf hinweisen, daB es, wie schon in der linearen Theorie, sowohl Frage- 
stellungen gibt, die sich prinzipiell nicht allein mit der Methode des Verschie- 
bungsfeldes bewaltigen lassen, als auch Probleme, bei denen man mit der Me- 
thode der Spannungsfunktionen allein nicht auskommt. Zu ersteren gehoren alle 
Probleme mit raumlich verletzten Kompatibilitatsbedingungen, zu letzteren alle 
Probleme mit Massenkraften ohne Potential. 
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Weiter gibt es Fragestellungen, die beiden Methoden zuganglich sind, bei 
denen jedoch eine Methode deutliche Vorteile bietet. Diese Bemerkung gilt auch 
fiir dreidimensionale Probleme. Entgegen friiheren Vorstellungen kann man auch 
hier oft vorteilhaft mit Spannungsfunktionen arbeiten, seitdem das Summations- 
problem der Eigenspannungen seine einfache Lésung mit Hilfe der Spannungs- 
funktionen gefunden hat ([25], [6], S.55ff.). Leider ist das dreidimensionale 
Randwertproblem bisher noch nicht allgemein mit Spannungsfunktionen geldst 
worden, doch lassen die Erfolge bei der Behandlung des Summationsproblems 
auch hierfiir Fortschritte erhoffen. 

Allgemein kann man etwa sagen: Besteht bei einem bestimmten Problem 
ein Unterschied in bezug auf den zur Methode des Verschiebungsfeldes und der 
Spannungsfunktionen gehérigen Arbeitsaufwand, so verstarkt sich dieser Unter- 
schied, wenn man zur quadratischen Naherung iibergeht. Die Vorteile der ein- 
zelnen Methoden bleiben beim Ubergang zu der héheren Naherung weitgehend 
erhalten. 

Beispielsweise behandelt man ebene Randwertprobleme oft vorteilhaft mit 
Hilfe der Airyschen Spannungsfunktion. Die Hauptvorteile sind 

a) Die Spannungsfunktion hat als Skalar die volle Symmetrie des betreffenden 
Problems und ist deshalb im allgemeinen von relativ einfacher Form (vgl. unsere 
Beispiele in §5 und §6). 

b) Man hat keine gekoppelten Differentialgleichungen oder Randbedingungen 
wie im Fall des Verschiebungsfeldes. 

c) Das Problem 1aBt sich sehr tibersichtlich als die Grundaufgabe der Bipoten- 
tialtheorie formulieren. 

Beispielsweise liege das ebene Problem eines Kérpers vor, der durch Rand- 
krafte und Massenkrafte mit Potential U beansprucht sei. Die Gleichgewichts- 


bedingungen lauten dann 
0; (0;; — U 0; ;) = 0: (89) 


Der Spannungsfunktionsansatz (@ = Massendichte) 
64=F,+eU, o,=F,+eU, oy=—F, (90) 
fiihrt dann in bekannter Weise auf die Gleichung ([3Z], S. 191) 
A AF = {(U) = — +—*" 9 AU, (91) 


wo f mit U als gegebene Funktion anzusehen ist. Die Randbedingungen lassen 
sich bekanntlich in der Form 
Pops), cron —1(s) (92) 


schreiben, wo g(s), h(s) als gegebene Funktionen auf dem Rand anzusehen sind. 
Die Gln. (91) und (92) definieren das Standardproblem der Bipotentialtheorie, 
das hinreichend ausfiihrlich im Schrifttum behandelt worden ist [7], [32]. Die 
Lésung dieses Problems liefert die Spannungsfunktion der linearen Naherung F°. 

Den Formeln des § 6 entnimmt man, daf jedes Problem, dessen lineare Lésung 
durch eine Airysche Spannungsfunktion beschrieben wird, auch in der quadra- 
tischen Naherung noch durch eine Airysche Funktion allein dargestellt werden 
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kann. (Das Entsprechende gilt nicht fiir die Spannungsfunktion der Torsion, 
vgl. § 5-) . 

Nach Aufsuchen der linearen Naherung laBt sich nun leicht zur quadratischen 
Naherung iibergehen. Wir schreiben 


A AF! = }' (F°) (93) 


und kénnen /’(F°) nunmehr als gegebene Funktion ansehen. Die jetzt zu erfiillen- 
den Randbedingungen sind offenbar 


Fi=0, 6F'/an=0. (94) 


Wir haben also, um die quadratische Naherung zu erhalten, noch einmal ein 
Standardproblem der Bipotentialtheorie zu lésen, wobei man sich weitgehend 
der Zwischenergebnisse der linearen Rechnung bedienen kann (Greensche Funk- 
tion, Matrixelemente bei Reihenentwicklung von F usw.). 


Man kénnte die Aufgabe auch mit Hilfe des Verschiebungsfeldes angreifen 
und wiirde auch hier zweimal ein Grundproblem zu lésen haben. Diese Bemer- 
kung zeigt deutlich, da8 sich der Vorteil einer bestimmten Methode in der linearen 
Theorie beim Ubergang zur quadratischen Theorie verstarkt. 


Die in §5 und §6 mit Hilfe der Spannungsfunktionen gelésten Probleme der 
singularen Schrauben- und Stufenversetzung sind etwa gleichzeitig auch von 
SEEGER und Mann [4] und von SEEGER und Buck [4] mit Hilfe der Verschie- 
bungen behandelt worden. Der Vergleich der beiden Methoden zeigt insbesondere 
bei der Stufenversetzung klar, daB im Falle von Versetzungen die Spannungs- 
funktionen das geeignete Hilfsmittel darstellen. 


Wir hatten auf Grund der Gln. (65) und (79) ohne gréBeren Aufwand auch 
manche Probleme mit stetig verteilten Schraubenversetzungen «,,(%, vy) oder 
Stufenversetzungen «,,(%, V), &,,(%, vy) behandeln kénnen. Diese Probleme ent- 
ziehen sich aus den obengenannten Griinden einer Behandlung mit Hilfe der 
Verschiebungen. In diesem Zusammenhang erscheint es sehr wiinschenswert, 
fiir die Spannungsfunktionen eine einfache physikalische Deutung zu haben. 
Der Grundgedanke unserer Theorie war die Vorstellung des nicht-Euklidischen 
nattirlichen Zustands. Diese laBt eine von SCHAEFER [16] gegebene Deutung 
der Spannungsfunktionen als besonders passend empfinden: Die Spannungs- 
funktionen werden als diejenigen Reaktionen (im Sinne von LAGRANGE) ange- 
sehen, welche das Medium im Euklidischen Raum festhalten, also seine Ent- 
spannung in den natiirlichen Zustand verhindern. ; 


Es erscheint uns weiterhin fiir viele Probleme ein Vorzug der hier beschriebenen 
Methode zu sein, daB sie durchgehend mit den Koordinaten des deformierten 
Zustands arbeitet [33], [28]. Tatsdchlich wird in den weitaus meisten Problemen 
die Beanspruchung des Kérpers als Funktion im deformierten Zustand gegeben 
sein, und in diesen Fallen ist eine Beschreibung des physikalischen Sachverhalts 
in den Koordinaten des deformierten Zustands unbedingt erstrebenswert. Nicht 
nur die physikalische Interpretation der berechneten GroéBen (etwa Spannungen) 
wird so einfacher, sondern auch die ganze Rechnung. So konnten wir z.B. auf 
die Jakobische Matrix, die in der Murnaghanschen Theorie [1] eine so wichtige 
Rolle spielt, bei unseren Rechnungen vollkommen verzichten. 
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Interessante Méglichkeiten ergeben sich weiterhin fiir die Behandlung der 
Temperaturspannungen, magneto- und elektrostriktiven Spannungen sowie der 
Konzentrationsspannungen. Diesen Fallen ist gemeinsam, daB man als Ursache 
der Spannungen eine spannungsfreie, sog. quasiplastische (oder eingepragte) 
Deformation ¢}, ansehen kann, welche dadurch, daB sie im Euklidischen Raum 
nicht ungehindert vonstatten gehen kann, zu gleichzeitigen elastischen Deforma- 
tionen ¢,, und Spannungen o;, fiihrt. Ware der Zwang, im Euklidischen Raum 
zu verbleiben, nicht vorhanden, so wiirde der K6rper nach Anlegen etwa des 
Temperaturfeldes sich spannungsfrei in den nicht-Euklidischen Raum begeben. 
Der so erhaltene nicht-Euklidische Zustand des Mediums spielt fiir die weiteren 
Uberlegungen die Rolle des natiirlichen Zustands der friiheren Paragraphen. 

Bezeichnen wir den Abstand zweier Punkte, deren gegenseitige Lage im defor- 
mierten Zustand durch dx* gekennzeichnet ist, im undeformierten Medium (ge- 
meint ist, vor Anlegen des Temperaturfeldes etc.) mit ds,, denjenigen derselben 
Punkte im nicht-Euklidischen und im deformierten Zustand wie friiher durch 
ds, bzw. ds, die zugehérigen Metriktensoren mit 0,,, g,;, 4,;, so ist offenbar der 
Eulersche Tensor ¢;) der quasiplastischen Deformation durch 


dst— dst =2e) dx" dx’, (95) 
derjenige der elastischen Deformation ¢,, durch 

is Ase eda da (96) 
zu definieren. Addition der beiden Gleichungen gibt 


ds*—dsi=2&,dx* dx (97) 
mit der Gesamtdeformation 
ef) = eh, + Ep. (98) 


Da diese Euklidisch ist, verschwindet der zugehérige Riemannsche Kriimmungs- 
tensor, dies liefert die Grundgleichungen fiir die Spannungsbestimmung zu 


nie te FIV, Date ae gD egl nip] tii) =), (99) 
wo jetzt jedoch 
Duin = — 2 (Grin + Emi) (100) 


zu setzen ist. «7, folgt im allgemeinen leicht aus den physikalischen Gegeben- 
heiten des Problems, z.B. ist bei gegebenem Temperaturfeld und einem Aus- 
dehnungsgesetz 


ele (101) 


mit %, V als Anfangs- bzw. Endvolumen eines Massenelements 
e41= (f(T) —4$f(T)?] Ont: (102) 


Man erhalt dieses Ergebnis leicht, wenn man mit den einer réumlich homogenen 

Ausdehnung (101) entsprechenden Verschiebungen uk —x,f(T) in die (nicht- 

linearen) Gleichungen fiir den Zusammenhang zwischen Verschiebungen und 
eee 
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Deformationen ((2], 5, 299) 
2ef = O,up + Ome — (0,u2) (0,47) (103) 


eingeht. : 
Offenbar ist es sinnvoll, nach Vergleich mit Gl. (23) —2éh12=imnx Oder, 
nach Gl. (34), weil Ap_.), in Gl. (28) nichts beitragt, 


eV, of = — a”, (104) 


als eine ,, Quasiversetzungsdichte‘ und diese als Ursprung der Spannungen an- 
zusehen. Diese Spannungen sind dann aus der Versetzungsdichte (103) wie die 
Eigenspannungen einer gewohnlichen Versetzungsdichte zu bestimmen. Bemer- 
kenswert ist, daB sich die Quasiversetzungsdichte im allgemeinen sofort als lokal 
ergibt, weshalb hier die Methode der Spannungsfunktionen derjenigen des Ver- 
schiebungsfeldes, die meist in den Anfangskoordinaten arbeitet, vorzuziehen ist. 


Der Begriff der Quasiversetzungen ist sehr niitzlich z.B. bei vielen Problemen 
der Metallkunde (([6], S. 32ff.). Hier entstehen beim Abschrecken von Werkstoff- 
proben starke Temperaturspannungen. Die in dem Medium vorhandenen (und 
eventuell neu zu bildenden) beweglichen Gitterversetzungen werden sich dann in 
jedem Moment so anzuordnen suchen, daB sie die Quasiversetzungen der Tem- 
peraturverteilung méglichst weitgehend kompensieren. Dabei geben die Gitter- 
versetzungen einer Schubspannung bei nicht zu hohen Temperaturen viel leichter 
nach (durch Gleiten) als einer hydrostatischen Spannung (durch Klettern). Dies 
ist der Grund dafiir, daB man keine sehr starken makroskopisch variierenden 
Schubspannungen (sog. Spannungen 1. Art) in einem gewoéhnlichen Festkérper 
anbringen kann (weder als Lastspannungen noch als Eigenspannungen); diese 
gleichen sich durch Nachgeben der Versetzungen, d.h. durch plastische Ver- 
formung, weitgehend aus. Fir Schubspannungen wird daher die Elastizitats- 
theorie 2. Ordnung hauptsachlich im Falle von mikroskopischen Spannungen 
(Spannungen 3. Art) bendtigt werden, also insbesondere fiir die Spannungen in 
der Umgebung einzelner Gitterfehler. 


Dagegen kénnen sich starke makroskopische hydrostatische Spannungen bei 
nicht zu hohen Temperaturen sehr leicht ausbilden, sie sind besonders einfach 
durch Temperatur- und Konzentrationsschwankungen zu realisieren. Sie kénnen 
so weit gesteigert werden, daB RiBbildung eintritt. Der Nutzen der nicht-linearen 
Elastizitatstheorie bei der Behandlung solcher Probleme liegt auf der Hand. 


Die Verfasser freuen sich, die vorliegende Arbeit ihrem verehrten Lehrer, Herrn 
Professor R. GRAMMEL, dem sie die erste Einfiihrung in die Mechanik und angewandte 
Mathematik verdanken, zu seinem siebzigsten Geburtstag widmen zu kénnen. Ferner 
méchten sie Herrn Professor E. Furs und Herrn Professor U. DEHLINGER fiir ihr 
Interesse herzlich danken, 


Anhang 
a) Wir stellen hier den Zusammenhang zwischen den von MURNAGHAN in [28 | 
beniitzten Elastizitatskonstanten dritter Ordnung /, m, 2 und den von uns be- 
nutzten Konstanten L, M,N her. Bei der ganzen Rechnung lassen wir Terme, 
die in o;; oder ¢,, von héherem als zweitem Grad sind, konsequent weg. Wir 
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beniitzen etwa die aus den Gln. (54) leicht abzuleitenden Beziehungen 


ey = 1, 0, + (3L +2M) 02? + (N— 6M) on 


3K 
4 9OKk*2— 4G? 
= 7G (on— 97K? oj) 
1 
Eyy — yal Oxy == 2M oy; Oxy aN \Gick Oya D2 Oxy) g 
wo 
ee 2G(1+7) 
3(1— 27) 


der Kompressionsmodul ist. Damit gehen wir in die nach MURNAGHAN folgenden 
Gleichungen 
oy =3K e+ (314 2m) ef + (n — 6m) ey 
Opp 2G Boy 2 Sj Egy tn WE ety, — Exe Ezy) 


/ 


ein (fiir ¢,,, ¢,, und ¢,, setzen wir 0,,/2G, o,,/2G und (611-45 1)/26, da 
diese GroBen nur in quadratischen Termen vorkommen). So erhalten wir 

OL +N =— (91+ n)/27K3 

6M — N =— (6m — n)/12K G? 


N =— n/8G°. 


Die Konstanten /, m, m verbinden, ebenso wie die L, M, N, die in den Endkoordi- 
naten genommenen (d.h. Eulerschen) Spannungen und Deformationen. Wegen 
des Zusammenhangs der /, m,n, mit den in MURNAGHANS Buch [1] beniitzten, 
sich auf die Lagrangeschen Spannungen und Deformationen beziehenden Kon- 
stanten, die von MURNAGHAN ebenfalls mit /, m, m bezeichnet wurden, vgl. SEEGER 
und MANN [4]. 

b) Wir zeigen, daB der Anteil mit «,, und «,, auf der rechten Seite von 
Gl. (79’") zu F1 im Falle der singularen Versetzung «,,—b06(0),%,y=0 nichts 
beitragt. Wegen F°~ylno ist AF°~y/o? und AFP~(x?— y?)/o4= (cos 2¢)/0?. 
Eine Partikularlésung der Gleichung 


A AF? ~ (cos 2¢) 6(0)/0” 
ist 
22 foe) 
F1 (0,9) ~fdy'cos29’' fdo'rd(o')lo’, P=e?+e0'?—2ee'cos(yp+ ’), 
0 0 


integriert tiber den ganzen Wertebereich von 9’ und g’. Die Integration tiber 9’ 
ergibt 


Vo?-+0?—20 0’ cos(p+ ¥) 
Q o'=0 


2a 
F4(0, 9) ~ [ dg'(cos 29’) 
0 


Dieses Integral verschwindet aber offensichtlich fiir simtliche Werte von 9 und 9. 
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A Note on the Exstencie of Periodic Solutions 
of the Navier-Stokes Equations 


JAMES SERRIN 


Communicated by C. TRUESDELL 


In an earlier paper [7] I established the stability of solutions of the Navier- 
Stokes equations provided that a certain Reynolds number was small enough. 
This result enables one to deduce from certain plausible hypotheses the existence 
of stable periodic solutions of the Navier-Stokes equations. In particular we 
shall prove the following result: 


LetW =¥ (t) be a bounded region in space, and let a flow velocity be prescribed 
at each point of the boundary of V. Assume furthermore that both V and the assigned 
velocities depend periodically on the time t. Then, under conditions 1° and 2° 
stated below, there exists a unique, stable, periodic solution of the Navier-Stokes 
equations in V which takes on the prescribed on the boundary of V. 


In view of condition 2° (see below) it is necessary that the boundary condi- 
tions be compatible with a flow of Reynolds number less than 5.7. The theorem 
may therefore be paraphrased by the statement that corresponding to suf- 
ficiently low assigned periodic velocities there exists a periodic flow to which 
every other motion eventually subsides. A simple example would be a fluid 
enclosed in a fixed container and stirred by a low speed bladed rotor. 

Another case of interest occurs when the assigned conditions are steady. Then 
the theorem asserts the existence of a unique, stable, time-independent solution 
of the Navier-Stokes equations taking on the prescribed velocities on the boundary 
of Y. In several long papers (e.g. [2]—[5]) the existence of steady flows has been 
proved independently of conditions 1° and 2°, though in these cases a con- 
siderable degree of smoothness has been required of the boundary data. Needless 
to say, the flows thus obtained may be neither stable nor unique, unless their 
Reynolds numbers are low enough. 


We now state conditions 1° and 2° envisaged in the above theorem and 
remarks, 


1°. To every continuous initial distribution of velocities over V there corresponds 


a solution of the Navier-Stokes equations satisfying the prescribed boundary condi- 
tions. 


2°. There is one solution whose Reynolds number Re=Vdly is less than 5.7. 


(Here V is the maximum speed of the flow during the whole time interval 0<t< 00, 


d ts the maximum diameter of V, and y is the kinematic viscosity.) This solution 
as equicontinuous in x =(x, y, 2) for all t. 
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Before proceeding to the proof of the theorem there are several remarks to 
be made. First, it is understood that the flows guaranteed by Condition 1° 
should be valid for all t20. Though this is mathematically a very stringent 
requirement, it is nevertheless quite plausible on the grounds of physical intuition 
and our knowledge of the behavior of other parabolic systems; moreover, from 
the proof it will be clear that 1° need hold only for initial velocity distributions 
whose Reynolds numbers are less than 5.7. Second, the number 5.7 in Condi- 
tion 2° may possibly be improved, as will be apparent from the proof. Finally, 
the author wishes to make it clear that he does not consider the above theorem 
to be equivalent to the usual kind of mathematical existence theorem: it is agreed 
that for certain types of boundary behavior 1° and 2° might not hold, and that 
(correspondingly) the limits of application of the theorem are not well-defined. 
On the other hand, in view of the plausibility of 1° and 2° there is every reason 
to expect that they will be proved in the future, subject to sufficiently smooth 
boundary conditions, and once such proofs are given the existence of periodic 
solutions indisputably follows. 

Now let v =v (a, t) be the velocity vector of the flow guaranteed by Condi- 
tion 2°, and suppose the assigned boundary conditions have period 1. We con- 
sider the sequence of vector fields 


DP) O(a, Ol. 5 een de oer. ac 


This sequence is bounded and equicontinuous by hypothesis, hence by the theorem 
of Arzela contains a subsequence which converges uniformly to a vector p(x). 
We assert that actually the whole sequence converges to @ (x). 

Indeed, if this were not the case we could find another subsequence converging 
to a vector (x) different from ~(x). But this is impossible, as we shall now 
show. For positive integers m and n, (m>n), set 

v(x, t) =v(x,i+m—n), t= 0. 


Then obviously v’ is a solution of the Navier-Stokes equations, and, according 
to the periodicity assumption, it satisfies the assigned boundary conditions. 
Now for any vector f=/f(#, t), let #(f) be defined by 


Hf) = 4% J fe, A)? dv. 
v(t) 
Then by Theorem 1 of [1] we have 
K(v' —v) Ss Hye", (1) 
where .% denotes the value of #(v’—v) at t=0, and 
é=— - (a v2/d? — V*) > 0 
since Re[v]<5.7. Because v’ as well as v satisfies ReS5.7, it is seen that 


Agate (QV )8 dv 2d°\V 2-—= Const. 
Ve 
KH (Pm — Pn) S Const. e~*”, (2) 


1 Similar theorems are already known, cf. the papers of LERay and KIsELEv & 
LADYSHENSKAYA cited in the references. 


Setting t= in (1) yields 
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and it follows that lim #(@~,,—©,) =0 (3) 


m, N—> CO 
(the region of integration in (2) and (3) is, of course, ¥(m) =¥V,). Now letting 
m and n tend to infinity through sequences of integers such that @,,— v,p,- Pp 
gives an immediate contradiction. This proves that lim @,, (x) =¢ (x). 
By condition 1° there exists a flow v*=v*(a, t) such that v* (x, 0) = (2). 
We assert that v* is a periodic solution of the Navier-Stokes equation. For set 
v(x, t) =v(a,t-+mn). Then we have, in the same way as inequality (1), 


Ht (v* — 0") < Hye", (4) 


where ¢ has the same meaning as before, and % is the value of #(v*—v’’) 
when t=0, that is #;=#(p—@,). Putting ¢=1 in (4) yields 


H[v*(@, 1) — Pnsa(#)] <4 — Pn) 
and letting 7—> oo gives, finally, 
HK |v*(a, 1) —ep(x)] =0. 
It follows that w*(a, 1) =p (x) =v* (x, 0), that is, v* is periodic. 


To complete the proof of the theorem it is enough to show that Re[w*|<5.7, 
for then by [1] v* is stable and unique”. But we have 


a (U* —v)=>0 “as 7—-co} (5) 


by virtue of Theorem1 of [1]. Since both wv and v* are equicontinuous, (5) 
implies (v*—v)—0 as too, and therefore 


Max v*(a#, 1) SV + 0/(2). 
Be 
v* being periodic, this in turn proves V*<V and Re[v*]<5.7. 


Note: This research was supported in part by the United States Air Force Office 
of Scientific Research under Contract AF 49 (638) — 262. 
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» That v* is stable follows from Theorem 1; that it is unique can be shown by the 
same argument used to prove Theorem 2. 
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The Evaluation of Collision Integrals, 
using Grads Representation of the Distribution Function 


E. IKENBERRY 


Communicated by C. TRUESDELL 


I, Introduction 
In the paper in which Grap [3] developed his thirteen moment approximation 
to the velocity distribution function for a rarefied gas he utilized certain Hermite 
polynomials whose properties he developed in a companion paper [4]. The 
Hermite polynomials, which are components of tensors defined by 


H (ae) = (—1)" eb" "4%, (1.1) 


are of total degree m in the three components %,, 2, 3 of a. The first few poly- 
nomials, as given by GRAD [4, eq. (14)], are 


HO —4, 


HY Ss 
TEND ze pe oe 
aj 44 ap? 

AY, = Hj Hj Hy — (%5 Oj» + %; Ong + %% O;,) - 
Since only these are needed in his thirteen moment approximation, GRAD com- 
puted the collision integral, which appears as the right hand member of the 
Maxwell-Boltzmann equation in its customary form, only for the above, obtaining 
by a lengthy development results exact for Maxwellian molecules and correct 
to third order terms for arbitrary molecular models. GRAD makes extensive use 
of a tensor notation which is not required in the present paper. Our definitions 
of the Hermite polynomials are equivalent to GRAD’s, although our work on 
the collision integral (Sections V and VI) is in the notation used by IKENBERRY & 
TRUESDELL [7]. 

Products of 3N one-dimensional Hermite polynomials were used by IKEN- 
BERRY [5,6] in a series representation of the distribution function in phase 
space. 

In this paper we wish to show how exact expressions for the collision integrals 
for any spherically symmetrical model may be obtained. We effect all integra- 
tions except those over the relative speed 2v and the scattering angle ®. For 
inverse power molecules, the integration over the relative speed will be effected. 
For rigid elastic spheres, all integrations will be effected. 
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II. The Hermite Polynomials 
Grap’s Hermite polynomials of degree » in three variables © =(%1, %2, %3) 
may be defined by the polynomial identity 


: : i 
(—A)* ct (eV) %e- 2" =n! yee = H™(r, x), (2.4) 
: ! 
where 
taxi t+ w+ 43, 
(4) (a) fa) 
i VS lea ries ee 3 Ox. ? 


rt = pis pir y's or 
Steg at acer eer Yo 


1 Stal te! tan, 
and the summation is over all partitions i=(7,, 72, 73) of m into non-negative 


integers. From (2.1) and the expression 


ns 


De ea 


it is readily seen that the three-dimensional polynomials H;(a#) are actually 
products of three one-dimensional Hermite polynomials H;,(%;): 


Hy (x) = H;,(%) Hi,(%2) H,,(%s) (2.4) 


vs 


where 
Fee VEN Ta pe y 
H,,(%;) = (—1)%@ Gael? ba (2.5) 
However, GRaAb’s contracted “Hermite” polynomials [3, eq. (5.6)| of degree 
m>1 are not products of one-dimensional polynomials. Following GrRapD [3, 
p. 363], we shall consider any given contracted polynomial as a certain linear 
combination of polynomials of the form (2.4). 


From (2.1) any number of the Hr, a) are readily obtained by noting that 


(7 -V) e?* =— (r-a) e3* 


> 


2.6 
(r-V) (rs x)= 73; aa) 
In order to obtain the few H;(a) given in (1.2) we need only 
De AAG iat) fa Ie 
HO(r, x) = (rx) =, ee 
H°2) bess aay (2.7) 
(Pe, a) —- (r : x) a ha — Uy V5 (X4 Mp 0; ;) ) 


H® (ya) = (r-x)®— 37 (r- a) = 151; 1 5 (Xi Xj Xp, — 3%; 0; x) - 


In the right hand members of (2.7) the summation convention has been used. 
ae symmetrizing the parenthesized expressions we obtain GRAD’s polynomials 
452): 

Integral and recurrence relations between the polynomials are most easily 
derived by means of their generating function. To obtain this, we divide each 
member of (2.1) by ! and sum the resulting expressions over all non-negative 
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values of x. Quite readily we obtain the generating function 


=p FEL. = Hi;(ax) rt Hn) : 
f(r, 2) =e74 eae, ur aw) (2.8) 


The easily derived definite integral, 
" e** tir, x) f(s, x) dx = (2n)t e's, (2.9) 
leads directly to the orthogonality and normalization relations for the H)(r, «) : 
Sed HO (x, a) HO (8, a) dar = (271)8 1! (1 8)" 3p, (2.10) 


Expanding the right hand member of this by the multinomial theorem (see (2.3)), 
and expressing the H™(r, x) and H™)(s, x) in the left hand member in terms 
of the H;(a) by means of (2.1), we readily obtain the orthogonality and normali- 
zation relations for the H;(x os 


foun x) H, (a) dx = (22)3 i! 0; ; (2.11) 


where Oi; — Oi jy Onis Shs 


III. Representations of the Distribution Function 


It will be assumed that, at least for some values of a non-negative para- 
meter h=h(a#, t) there exists a uniformly convergent series representation [3, 


eq. (4.6) | 


he _gaace wa 0; (a, 2) 
i ead NESE BAUME <i 
of the distribution function F=F(a#,&,t). Here & is the velocity of a molecule 
relative to some standard frame of reference that is the same for all molecules, 
and 


c=——C (3.2) 


where C = C(a, t) is the velocity of a selected coordinate system for the mole- 
cule at a. As a consequence of (2.11) and the assumed uniform convergence of 
the series in (3.1), we find readily that 


b; = b;(x, t) = J FH; (he) de. (3.3) 


The 0; divided by the number density ” give GRAD’S a;. GRAD [3, eq. (4.11) | 
gives expressions, in terms of moments of F about the mean velocity , for those 
a, for which 7,+7.+7,=0, 1, 2, 3, and 4, when C = € and when h?=m/kT, where 
m is the mass of a molecule, k is the Boltzmann constant and 7 is the tempera- 
ture. It is readily verified, as stated by GRAD [3, eq. (4.12) ], that taking C,=€, 
implies a;=0, and that taking h?=m/kT implies a\)—0. The calculations of 
the present paper are made with no specification regarding C and h?. 


IV. The collision integral; Notation 


Solution of the general initial value problem for the velocity distribution 
function F by the transport equation method involves calculation of the collision 
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integral oy Ieee 
C(Q) =fdcfdet [de fvbdbdQrr* (4.1) 
(oe) [o-e) 0 0 


for various polynomials Q =Q(a, €, f). Of the permissible alternative forms for 
dQ we select s0=0" + 0'— 0* 0, (4.2) 


which shows the invariance of C(Q) when e*<e and c*’<c’. 

In the above and in the following developments we adhere closely to rather 
standard notations, in particular to the notation used by IKENBERRY & TRUES- 
DELL [7], from whom we quote the following paragraph [7, p. 15] for the con- 
venience of readers of the present paper: 

‘Quantities associated with the second of a pair of molecules are denoted 
by an asterisk; quantities associated with the outcome of an encounter, by primes; 
and the accents when attached to function symbols refer to their arguments: 
F'= F(&’), etc. We set 


20 = Eb = CF 6, 2a Co md) (4.3) 


Conservation of momentum and energy requires w=w’, v=v’. Let e be the 
angle between the plane of » and v’ and the plane containing v and a direction 
fixed in space. Suppose § fixed; then all possible encounters will be included 
by letting ¢ increase from 0 to 27; in this process v’ rotates around w with fixed 
magnitude v’=v while subtending a constant angle ® from v. To include all 
possible encounters experienced by a molecule of velocity § at x, we integrate 
over a plane normal to vw and passing through a. In this plane, let b, ¢ be polar 
coordinates. Then 0 is the distance from & to the initial asymptote of the second 
molecular trajectory. A knowledge of the law of intermolecular force enables 
us to calculate 6 as a function of ® and v.” 

For a spherically asymmetric molecular model, } is a function of ¢ as well 
as of ® and of v. Since we are considering only spherically symmetric molecular 
models, as implied by the above paragraph, we shall write 


db| __ 


EON ree 


B(v, ®) =v 1 (v, B) (4.4) 


where I(v,®) Csc@ is the differential scattering cross section. 

In the next section we perform some of the integrations in the evaluation 
of the collision integrals C(Q) for selected Q. In Section VI we show how the 
integrations for rigid elastic spheres and for inverse power molecules can be more 
readily effected. ; 

V. Evaluation of the Collision Integrals 


Substituting the series representation (3.1) into the expression (4.1) for C(Q) 
and effecting a simple change of scale by replacing he by ¢, h?de by de, hv by 2, 
etc., we obtain 


C(Q) = 2 5b; Ri; (Q) (5.1) 
where J 
Rij(Q) = gis [de [der [ae [aw Bp Pert cig eae ee) 
‘ P " uz: 


ce) 0 0 
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We take Q=H;, (he). Designating by R,;(k) the value of R,j(H;,(e)), we 
see by means of (2.8) and (5.2) that R,;(k) died by k! is on coefficient of 
p'qir® in the power series representation of 


nena om [te fact fee fan o(n A (Atom) 5 


MCP | ere er & Cee 2 


(5.3) 


Here 
LT Pa Te (5.4) 


For brevity, we write R in place of R(p, q, r). 


The quantities b;, defined by (3.3), are expressible in terms of moments of 
F and may be regarded as representing gross properties of the gas flow. F.g., 
those of lowest index are certain combinations of the density, the mean velocity, 
the temperature, the pressure tensor, and the flux of energy vector. By (5.2), 
quantities R;;(Q) are determinate from the intermolecular force law, without 
knowledge of F’. To find the R;;(k), it suffices to calculate formally the power 
series representation of R. Therefore, explicit evaluation of all collision integrals, 
for any spherically symmetric molecular model, 1s equivalent to conversion of the 
definite integral (5.3) for R into a power series. To this task we now address our- 
selves. 

We start by transforming to w and v in place of € and e* as variables of 
integration in (5.3). By means of (4.3) and the relations dedce*=8 dw dv, 
w=wu’ and v =v’, we obtain 


k= MS" fae fae fue faon(y ®)e en vnwt x 


x 08 BREW (Cosh (I. w') — Cosh (r- B)}, a 
where 
S=q—-PD, 
pe: i 
Since [dw ew +8810 mh et, (5.7) 
(5.5) may be written in the form 
Late fooalb eon 
bes 5 0 (5.8) 


x {Cosh (r- v’) — Cosh(r-v)}. 


We next perform the integration over de. Taking for ¢ the angle between 
the planes determined by v and wv’ and by v and r, we have, from a cosine law 
in spherical trigonometry, 


r-v'=r-vCos® + rv Siny Sin® Cos «, (5.9) 
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where y=(r,v) and ®=(v,v’). Then, by means of the definite integral 
22 
{deer — 2a TitA), (5.10) 
0 


where I,(A) is the modified Bessel function of order zero, we obtain, for any 
spherically symmetric molecular model, 


Aye p= 7? eo \ -v+ts-0 
Raid fav [a@B(y, Pe x 
eS 9 (5.11) 
x {Cosh (1 - v Cos ®) I(A) — Cosh (r- v) 5, 


ere A=rv Sinz Sin®. (5-42) 


We proceed by representing the integrand of (5.11) as a series in the com- 
ponents of p,q and r. First we write 


a A\O gh 8 a8 
ak Ege =) ae = (5.13) 
x 
ner ba ok aera (5.14) 


the sextuple summation being on all triplets @ = (%, %, #3), B = (1, B2, Bs) of non- 
negative integers. We also write 


2¥ w¥ Cos?’ D 


Cosh (r- v Cos®) = 2 = ay! (5.15) 
Since co 42h 
RO = ss 5.16 
Iii ervas (5.16) 
and 92 = {72 v2 — (r-v)"} Sin?D, (5.17) 
we have, by using the binomial expansion 
= ra ' (—1)° x9 y° 
(x—y) MO a, (5.18) 
the expansion 
I, (rv Sin x Sin ®) = y (=A)? (sv)? oP? (in @yzetee (5.19) 
oeno 4°** dl el (d+.e)! 


The implied changes in the order of summations here, as in later expansions, 
are conveniently legitimate without consideration of convergence since we are 
here interested only in obtaining the coefficients in the series representation of R, 
not its evaluation for numerical values of the components of p, q and r. 


By means of (5.19) and the trinomial expansions 
t=O) 5 


e2¢€ 5.20 
(7° 10) 2S a Dahil ined ( ) 


* Only for the triplets a, B, y. §, €, A, Uv, e, T,X, and w, identified by Greek letters, 
do we use symbols a= o,+a ,+0a,, etc. 
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we obtain 


€ $6+2€ 26 2e 
I, (vv Sin y Sin®) = (=1)* (22)1r* v©  (Sin@)22+28. 
a 2 4°F®§1 (2e)! e! (+e)! ven) 23) 


In view of (5.15) and (5.21) we write the expansion of Cosh(r-v) in the form 


)o (2e) lat vt25t2e,, 26 perytze 


Cosh (r - v) ue és 
=) d+e,0 aaa WEEE (5.22) 
Finally, substituting (5.13), (5.15), (5.24) and (5.22) into (5.11), we obtain 
ie. -8 (—1)*** D&E Pe p* q® 2 (vt8+e) 
yes 
R=e > a!B! (2y) 8! (2e)! (5.23) 
where 
Dy cae d d®B en pRt2yt2€ y26 
vs LER i f (+: ®) = Hoy 
(5.24) 


xf Cos B)?” (Sin P)?9+?* — 6, ,, of. 


The summation in (5.23) is over the five sets of triplets a, B, y, 5, and € of non- 
negative integers. Letting 


K5(v) Sei d® B(v, D) {Cos’” D Sin? *D — 6, oh, (5.25) 
we have 
ae = eels ee fea Gur OO ATT 6 720) (5.26) 


Finally, using polar coordinates (v, #, gy) to integrate over the surface of a sphere 
of radius v, we obtain 


[oe) 

2e)! ae 

Doge - Tey Ansavene [dv Kyt*(2)e Papen Ain 4 one oe (5:27) 
0 


where 


2a a 


A, =4 | dp f 49 (Cos 8)” (Sin 3)*+%** (Cosp)” (Sing) 
0 0 


eee lag ec) 


(5.28) 


We are able to effect the integration over d@ in (5.25) only for rigid elastic 
spheres, and that over dv in (5.27) only for rigid elastic spheres and for inverse 
power molecules. More direct methods of performing the integrations for these 
molecular models are given in Section VI. 

We now complete the expansion of FR as a series in the components of p, q 
and r. From (5.4) and (5.6) we express #?— JT? in terms of p,q and r. Then, 
using expansions of the form (5.13), we obtain from (5.23) 

AeA? 


«+B 
hoe Yee ee epere ca I (5.29) 
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where o = (a, B, 27,5, 2€, A, u,v, QT, xX); 

o! =a! B! (2y)!8! (2e)!Alply! pe! ex! 
i=a+2A4+T+YX, 
j=P+2e+e7X, (5.34) 
k=2y+26+2e€+2v+ 9+: 


(5.30) 


and where 


The o under the summation sign indicates a 33-fold summation over non-negative 
integers. In (5.34), addition and equality of triplets are defined as in standard 


vector notation. ey 
Since R,;(k) is equal to k! times the coefficient of p'qJr* in the expansion 
of R=R(p, q, Yr) as a series in the components of p, q and r, we see immediately 


from (5.29) and (5.31) that 
fe Be ie” 


= a+B 
R,;(k) =k! 2 ith petty) DYES (5.32) 


where the summation is over all o satisfying (5.31), i,j and k being fixed. 


The terms in the right hand member of (5.32) may be symmetrized in @ 
and £ by replacing (— 1)* by $(— 1)*+ 3( 1)®. This step, which may be justified 
by examining (5.28), may be traced back to the symmetry in € and e* of the 
integrand in the general collision integral (4.1). The vanishing of A,, for many 


w greatly reduces the number of non-vanishing DE (5.52). 


We summarize our general result. Six of the eight quadratures needed to 
determine all collision integrals for any spherically symmetric model have been 
effected explicitly. By (5.32), Ry;(k) is expressed as a finite sum of the quantities 


Dae which are given as double integrals by (5.26) and (5.27). 


VI. Special Molecular Models 
A. Rigid Elastic Spheres. For rigid elastic spheres of diameter d, for which 
[2, p. 168; 3, pp. 342—343] 
b=dCost@, 


= Origin (aioe engee (6.1) 
Biv, B) =vb 75 |= ee Sin®, 


the integrations over d® and dv in (5.25) and (5.27), giving K%(v) and pie are 
easily performed. A better procedure, however, is to observe that de d@® Sin ® 
is the differential element of solid angle in the direction of wv’, in a coordinate 
system with v as polar axis. Taking polar coordinates (v’, y’, e’) for v’ in a 
coordinate system with r as polar axis, we readily find that 

20 1% 22 na 

J de fd®Sin®Cosh(r-v') =fde' fdyz' Sinz’ Cosh (ro Cos x’) 

0 0 0 0 

Sinh vv 
0 
VU 


(6.2) 
=4 


Hence we obtain, from (5.8), (6.1) and (6.2) 


a 2a% 2» 7: epee Sint 
R= 2e el ” fave + "9 Cosh(r- a). (6.3) 


‘a 
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An integration over the surface of a sphere of radius v may now be performed, 
but then a series expansion would be required before integration over dv. The 
final result is more readily obtained by expanding the integrand in (6.3) in terms 
of p, q and r and integrating termwise over dv. Using (5.13), (5.15) and (5.20), 
we obtain 

Pa (4¢ pee 


Rah ATT peg p27 6. 
h Ay alplayr Pa Fy 
where 
= a Te eV ohvl a ere 
(2y+1)! ¥ (6.5) 


a? B (2y)! 7! 
Zi sae Fs +y +2) { our Ap— Ansar}. 
A,, having been defined by (5.30). 
We now express #?— J? in the exponent in (6.4) in terms of p, q and r, by 


means of (5.4) and (5.6), and expand the rsulting exponential expression in terms 
of the components of p, q and r. Using expressions of the form (5.13), we obtain 


AN a he ate eee wee 
a ya a Bee won pair (6.6) 


where, here, 
o = (4, 8, 27,0, u,v, 9,7, x), 


o! =a! B!(2y)!Alplvie!t!x!, Sie 


and 
t=a+2A4+7+4+Y, 
j=B+2p+e+x, (6.8) 
k=2y -2Y +p --t- 
The o under the summation sign in (6.6) indicates a 27-fold summation over non- 
negative integers. 
Finally, from (6.6) and (6.8), we see that, for rigid elastic spheres, 


Atputr 
ay eee (—1)** : a+p 
R,;(k) =k! 2 ae ae pe (6.9) 


where the sum is on all o satisfying (6.8), i,j and k being fixed. In (6.9), we 
may symmetrize in a and B by replacing (— 1)* by 3(— 4)t 48 (ar. 

On comparing (6.6) to (6.9) with (5.29) to (5.32), we see that the method of 
integration for elastic spheres is equivalent to effecting summations over 6 and 
€ in (5.29). 

B. Inverse Power Molecules. For a repulsive force law varying inversely as. 
the s** power of the distance d between molecules, it may be shown that [J, 
pp. 153—157; 2, pp. 170—171; 3, p. 342; 8, pp. 40—43] 


B(v, ©) =v" f, (®) (6.10) 
where 
oe (6.11) 


O= 
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Quite evidently, the special form of (6.10) permits us to integrate over dv 
in (5.24). However, the remaining integration, over d@®, can be effected only 
by lengthy numerical calculations. 

Substituting (6.10) into (5.24), and integrating over dv, we find that 


p® Beaty totes PEP S) Apsayrne Sp *(m) (6-12) 


where 


Se(n) = gat [ EPh(®) {Cos?? @ Sin?’ O — 6, o}. (6.13) 
0 


For the Maxwellian gas, for which s=5, =O, there is no especial simpli- 
fication in these formulae. However, by means of arguments used by GRAD 
[3, pp. 360—361] we can deduce from (5.2) that, for Maxwellian molecules, 
Rj; (kt) =0 unless 11+ %)+%3+41+J2+)3='y + ko+hg. Thus, for Maxwellian mole- 
cules, the sum (5.1) giving C(Q) becomes a finite sum, for Q any polynomial. 
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Spherical Means and Radiation Conditions 


CALVIN H.WILcOx 


Communicated by A. ERDELYI 


Introduction 


Nearly fifty years have passed since SOMMERFELD introduced his radiation 
condition [9]. During this period it has become customary to use the condition 
in formulating and solving the boundary value problems associated with the 
diffraction of time-harmonic waves. The radiation condition is satisfactory from 
the mathematical viewpoint in that it leads to boundary value problems having 
unique solutions [7]. However, the physical reasons usually advanced for adopt- 
ing it, rather than some other condition, are far from convincing [9, 10, p. 190]. 
Our purpose here is to provide a more satisfying foundation for the use of Som- 
MERFELD’S condition by deriving it from other facts concerning wave propagation 
that are both mathematically demonstrable and evident to physical intuition. 


The diffraction problem for time-harmonic waves is usually introduced by 
considering wave sources having a time-harmonic, or sinusoidal, time dependence 
and assuming that the wave generated by these sources has the same time depend- 
ence. This assumption is attractive because it relieves us from the necessity of 
considering the initial values of the wave, thus replacing an initial-boundary 
value problem by a pure boundary value problem. However, we pay for the 
simplification in that the wave is not uniquely determined by the conditions 
imposed, on account of the possible existence of standing waves. To see why 
the uniqueness has been lost we must examine more carefully our assumption 
that the wave is time-harmonic. We can expect this to be strictly true only if 
the sources have acted for an infinite length of time so that any transient waves, 
produced when the sources began to act, have vanished. If we assume that the 
sources begin at a finite time then we must specify the initial condition of the 
wave and the wave is uniquely determined thereby. It is not in general time- 
harmonic, but physical intuition suggests it will differ from a time-harmonic 
wave only by a transient wave. Thus we might find the desired time-harmonic 
wave by solving the suggested initial-boundary value problem and passing to 
the limit of large time. Of course, if we can carry out this program then the 
problem is solved and we have no need for a radiation condition. Actually, this 
is possible in only a few cases. Hence we seek to characterize the time-harmonic 
wave directly, as the solution of a boundary value problem, without reference 
to the more difficult intitial-boundary value problem. This is accomplished here 
by showing that SOMMERFELD’s radiation condition for the time-harmonic wave 
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is a consequence of a simple and easily demonstrable property of the solution 
of the initial-boundary value problem; namely, that it represents a wave pro- 
pagating outward from the sources with unit velocity. 

A similar program for finding radiation conditions has been proposed by 
J. J. Sroxer [11, 12, p. 175]. However, in his work it is carried through only 
for two problems, those of surface-gravity waves in an infinite ocean and waves 
on a running stream, where closed analytical expressions for the solutions of the 
complete initial-boundary value problems are available. 

In § 4 the initial-boundary value problem for time-harmonic sources is formu- 
lated and the solution is shown to consist of a time-harmonic wave plus a transient. 
In § 2 the basic properties of the spherical means of wave functions are developed 
and applied to show that any time-harmonic wave function is the sum of a radia- 
tion function and a standing wave. §3 and §4 are devoted to the derivation 
of SOMMERFELD’s and several equivalent radiation conditions. In § 5 correspond- 
ing results are deduced for the vector diffraction problem. 


§ 1. The Initial-Boundary Value Problems 


Solutions u(p, ¢) of the wave equation 


_ Pulp, t) 
Au(p,t) = ag 
represent disturbances propagating with unit velocity. Indeed, it is to this 
fundamental property that the equation owes its name. The property finds an 
analytical expression in 


Theorem 1: Let u(p, t) be defined and satisfy the wave equation in the space- 
time domain D defined by 
PS TESE Uy, SOs 


where v 1s the distance from p to a fixed point O, and r,=0 is a constant. Moreover, 
let 
du(p, O 


u(p,0+) =0, pe at 2 [OP Vi te: 


Then u(p, t) vanishes identically in D. 
A less precise, but more suggestive, statement of the theorem is that if a 
disturbance is contained in the sphere ©(O, 7), with center at O and radius Tan 


at time ¢=0 (i.¢., w and 0u/dt vanish outside G(O, 79) for t=0) then at any later 
instant ¢ the disturbance is contained in © (O, %+2). 


The proof of Theorem 1 is effected by showing that any solution w(, 2) of 
the wave equation in 2 can be expressed in terms of its initial values 


“(,0+) =F), 2H(b-OF) _ Gi) 
by Potsson’s formula 
u (6,1) =4G(p,t) + eF(p, o}, 


the formula being valid at every point (£, t) in ®. Here F (p, t) is the mean value 
of F over ©(#,t). In particular, if F and G vanish then u(p, t) is necessarily zero 
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in Q. The proof that w(p, ¢) is given by Porsson’s formula can be given by a 
classical method due to LiouviLLE; see [13, pp. 212—213]. It is not reproduced 
here. 

Now let us consider the propagation of waves in an exterior domain & in 
Euclidean space; that is, a domain & that includes all points lying outside of 
some sphere. The boundary © of & is then a bounded set which we assume, for 
simplicity, consists of a finite number of smooth surfaces. The physical process 
that we wish to consider is the scattering, by objects represented by the boundary 
©, of the wave generated by time-harmonic sources of density 


(d, k) gah. 


which begin to act at time =0. The nature of the scattering objects is assumed 
to be described by a Dirichlet or Neumann boundary condition on ©. Thus, if 
we assume the initial state of the wave medium is one of rest, the process is 
described by a wave function (f, ¢) that satisfies 


Ap(p,t) -22@4 — Hp, ae, pin B, t>0, 


(1.1) 9(p,0-+) =0, op (bi04) =) pin &, 


OP, t) =O oF ee von t= 0; 


where 0q/0n is the usual normal derivative of m on S. It is also assumed that 
the sources are confined to a bounded portion of ¥%, so that /(f, k) vanishes 
outside a sufficiently large sphere. 

The wave function (f, ¢) is uniquely determined by conditions (1.1). To 
show this we make use of Theorem 1 and the energy theorem for the wave 
equation. The latter asserts that if u(,7¢) satisfies the wave equation in the 
spatial domain ¥ for a<t<b then 


(1.2) edn Ge er shie aemeh NEE Sig a a0; 


where © is the boundary of ¥. It is easily verified by differentiating under the 
integral sign on the left and using the wave equation and the divergence theorem. 
Now, to prove the uniqueness of y(f, t) let y’(#, ¢) be any other solution of (1.1) 
so that u(p, t) =@(p, t) —q'(p, t) solves (1.1) with {(p, k) =0. We complete the 
proof by showing that u(p, T)=0 in ¥ for any fixed T>0. First of all, if 
S(O, 7) is any sphere containing S and outside of which /(f, k) =0 then both 
y(p, T) and g’(p, T), and therefore u(p, T), are zero outside of S(O, 7+ T), 
by Theorem 1. Now apply (4.2) to w in the domain &7 consisting of those points 
of B that lie inside S(0,7),+ 7 +1). The surface integral in (1.2) falls into two 
parts, one over S and one over ©(O,7)+ 7 +1). The first vanishes because u 
or 0u/On is identically zero on S. The second vanishes because w is identically 
zero on G(O,7%+7 +1) foroS<i<T-+1. Thus 


ee asnaais C2 ri 


ot 
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On integrating from 0 to T and using the initial conditions we get 


24 (ou)? = fi 2 ay dV =0. 
fe a a (ape VN ( a t=0+ 
Bp Br 
Hence u(p, T) is a constant which must be zero because of the initial conditions. 
Thus u(p, T) vanishes in Bp and therefore in all of &. 

Our next objective is to show that the wave function g(f, ¢) can be written 
as a sum 
(1.3) p(b,t) =u(p,t) +e" uA), 
where u(f, ¢t) is a transient wave, 

jim u(p,t) =0, 

and u(p) is independent of ¢. It is clear that such a decomposition is unique if 
it exists. The existence of the decomposition can be shown by considering a 


related initial-boundary value problem for the domain ¥&, the pulse problem 
[3, 4,5]. It is defined by 


A Go (P, t) ee = 0,7 pn So (20, 


(1.4) poh ot)=0, Seal OF) _— _J(p,h), pin B, 


Pol, t) =O or Se oe OnTSs— i= 0. 


These conditions are formally equivalent to requiring that qo(p, ¢) be defined 
for all ¢, vanish identically for <0, and satisfy the differential equation 


2 
(1.5) Ayo(p,t) — 2 Pel?) — 4(p, 2) d(), p in B, allt, 
while the condition on © is unchanged. Thus qo(A, ¢) describes the wave due to 
impulsive forces acting at time ¢=0 in a medium that is initially at rest. It is 
this interpretation that suggests the name “pulse problem”’ for (1.4). The same 
considerations that were applied to problem (1.1) above show that the pulse 
wave function go(P, ¢) is uniquely determined by conditions (1.4) and that ,(f, 2) 
vanishes outside ©(O, 7)+4¢) for each fixed t= 0. 

The relationship between the pulse wave function and (f, ¢) can be obtained 
by the principle of DUHAMEL, or superposition principles [3, 4, 5]. The result is 


Theorem 2: 
pp, t) =I 90l6,t— ner dee epAnsateew. 
This equation is formally evident from (1.5). To verify it let us define y’(A, ¢) by 
(=f polb,t— eae 


and show that it satisfies (1.1). It must then coincide with y(p, t) by the uni- 
queness theorem. It is clear that the boundary condition satisfied by g on © 
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implies the same condition for g’. Next 


t 


t 
en = 9 (p, 0+) gone O9o(p, t—T) e-tkt dz =/ Se ae e-tkt dy 
0 


ot t 
0 


whence g’(f,0-+) and Hale vanish. Finally, 


coe ot 
t 


o? tis —tkt ae : 
tf PolP Poth ide = Fpl ble LA Gp, H), 


 p'(p, t) ab! OD (P, 0+) ew tht 


0 
which completes the verification. 


We can now discuss the decomposition of g(, t) into a transient and a time- 
harmonic wave. Notice that under a simple change of variable Theorem 2 
becomes 


p(p,t) = aoa 1) dt. 


Thus if m(p, t) has a decomposition (1.3) with w(, ¢) a transient then 


(1.6) J e'* o(p, t) dt 
0 
necessarily exists and equals (A): 
(1.7) u(p) ee: t) dt, 
and consequently 
(1.8) u(p,t) = —e-** f e** oo (p, t) dt = — J e'** po(p,t +8) ds. 
t 0 


Conversely, if the integral (1.6) exists then p(p, ¢) has a decomposition (1.3) 
where u(f, t) is given by (1.8) and hence is a transient. Moreover, if the last 
integral can be differentiated twice under the integral sign then (A, ¢) satisfies 
the wave equation 
pei 
0 


because gp satisfies this equation. It then follows from (1.3) and the equation 
for w that w(p) satisfies 


(1.9) Au(p) + k?u(p) =f(,%). 


Finally, it is clear from (1.8) that «(f,¢) and therefore u(f) must satisfy the 
boundary condition 
Ou 
id 


u(p) =O or EE oy; PODS: 

Our physical interpretation of the initial-boundary value problem suggests 
that (1.6) must converge and define a twice differentiable function. Indeed, if 
Q is the whole of space, so there are no scatterers, then @ (Pf, t) is given by 
Poisson’s formula. Hence, for each fixed p, go(p, 4) actually vanishes for all 
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sufficiently large ¢ and (1.6) is a proper integral. In the general case @p (p, #) 
consists of a direct wave, given by Porsson’s formula, plus a reflected wave or 
echo. Physical intuition suggests that the latter is due toa “trapped”’ wave 
that propagates along the scatterer © and radiates at each instant a wave with 
amplitude proportional to its instantaneous amplitude, which suggests that 
@o(P, t) decreases exponentially. However, present knowledge of the initial- 
boundary value problem is insufficient to prove this. . 

In the sequel we shall assume that (1.7) converges to a solution of (1.9). We 
shall call the problem of determining u(p) the diffraction problem for the scat- 
terers G and sources f(p, k). Our object is to characterize u(p) directly, without 
reference to q(p, 4). Now we know that w(f) satisfies 


Au+Ra=f m% 


(1.10) sep) San Wee ould) _ o, J ates 


but that these conditions fail to determine a unique function. The simplest 
example showing this is perhaps the function 


u'(p) = 


which satisfies the Helmholtz equation 


sin Ry 
Y 


Au’ + k?u’=0 


and vanishes on the spheres of radii wn/k about the origin, m =1, 2,3,.... On 
the other hand, RELLIcH [7] has shown that there is at most one function satis- 
fying (4.10) and SOMMERFELD’s radiation condition. Hence, we will have charac- 
terized u (f) if we can show that it must satisfy this condition. This is accomplished 
below. We begin by developing some properties of the spherical means of wave 
functions that are needed in the proof. 


§ 2. Spherical Means and Wave Functions 


If u(p) is a point function on a domain & in Euclidean space we denote by 
u(p,7) the mean value of u(f) over the sphere G(f,7) with center at p and 
radius 7; thus 

is 1 
Padi See 
(2.1) u(p, 7) oe udS. 
S(p,7) 

If u(p) is continuous in B then %(, 7) is defined for all pairs (p, 7) such that 
S(p, 7) is contained in B. For our discussion of the properties of these spherical 
means it is convenient to identify the points p with their position vectors # = fe) p 
relative to a fixed origin O. Then to each point function u(p) there corresponds 


a function of # which we denote by the same letter: u(p) =u(x). With this 
notation (2.1) becomes 


(2.2) ii(w, 7) = [uw +r) dQ, 


where {2 is the unit sphere, dQ is the element of area on Q, and # is a unit vector 
defining a unique point on Q, 


| 
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We shall need the following result, due to DarBoux. If u(«) is of class C? 
in a domain & then # (a, 7) is of class C? in its corresponding domain of definition 
and satisfies there DAaRBoUX’s equation 


where A is the Laplacian relative to a. For a proof see [2, p. 411]. Notice that 
DaARBOUX’s equation can also be written 
a? (7, it) 

or? 


(2.3) A(ru) = 


Let & be an exterior domain and let 7)(a%) be the least radius 7 for which 
©(#,7) contains the boundary of ¥. Then if w(a) is defined in B, %(a, 7) is 
defined for every # and 7>7,(a). Now let g(a, ¢) be a wave function, 7.e., 
we 0? p (a, t) 


Ay (x, t) er , 


for x in 8 and ¢>0. Then the spherical mean @ (a, t, 7) satisfies (2.3) for all a, 
7>Y7 (a), and ¢>0. Moreover, we can apply A under the integral sign in (2.2) 
and use the wave equation for g. The result is 


Pry) _ @(r@) 


ot? Cy ae 
719 (x,t,r) =F(r —t) + Gir +2), 


whence 


for each fixed #. The functions F and G can clearly change with x. Thus we 
may write 
(2.4) (a, t,r) = F(x, v—t)+ G(x, 7+?) ; 


if 


F and G can be determined from the relations 


PNGAAGT SDs D4 aed @ fetta 
F'(a,7—t) (+ m)re az | (& a) PP (e rr, t) dQ, 

2.5 @ 
oe eS WR ae ee et gael tare eno 
(7 +4) = (oot az)7% oa et TERETE Lake, 

Q 


which makes clear that F(a, t) and G(a, t) are functions of class C?. (We have 
written F’’(a,t) for oF (x, t).) It is also evident from these relations that 
F(a, t) and G(a, ¢) are wave functions for all # and for all ¢ and t>7, (ax), respec- 
tively. We shall apply these results to the pulse wave function gp, in the next 
section. 

Next consider a Helmholtz wave function u(x), 7.e., a solution of the Helm- 
holtz equation Aule) = ule) 20, 
defined and satisfying the differential equation for x in ¥. This equation, together 
with DARBOUX’s equation (2.3), implies that the means u% (a, 7) satisfy 

a? (vu) 
or 


+ h(rit) =0 
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for all a and r>7,(x). Hence 
vii (x, 7) = Ae **? + Bei’, 
where A and B are independent of 7 but may depend on a. We shall write 


% ; i evkr 
(2.6) ii (ee, 1) = — —— hy (@) — tea), 


the constants ee being introduced to simplify subsequent formulas. From 


(2.6) we have os 
(=. ee ik)ri =e (x), [ 


or \ OY 
whence 
ikr é * 
Uy (x) =< [ (Go— th) rue + 7%) dQ, 
(2.7) ae | 
te (ae) = ee [G, ik)ru(e + rr) dQ. 
Q 


By differentiating these expressions under the integral sign and using the Helm- 
holtz equation for u we find that (a) and u,(ax) solve the Helmholtz equation 
for all x; that is, they are entire Helmholtz wave functions. 


The functions uw, and uw, also occur naturally in connection with the problem 
of representing w in terms of layers of simple and double sources on the boundary 
S of B. To show this, let p be any point of ¥ and apply the divergence theorem 
in the form 

i na = pNP eas 
[ {ude vAuldV = ‘i {u aa v5 }ds 
D Boundary of D 
ikR 

to u(q) and v(q) = > where RF is the distance from # to q, in the domain D 
lying between © and a surface S’ containing & and # but excluding G(f, «) 
where « is small. The integral over D vanishes. Moreover, the integral over 
S(p, %) equals 4 u(p), as we can see by noticing that it is independent of « 
and tends to 42 u(p) for «+0. Thus we have 


HO) aa eee) oe 
S 


Af) Bos 


where ©’ is any surface containing p and ©. On taking S’= S(p, 7), so that 
R=r and n=r the last integral becomes 


A Ou |/ erkr 8 / etkr 
4% he | Y +) : al Y jas. 
S(p,7) 


Since 
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this integral equals 


(ety 2 [rue tri) aa) — es I ( q ik) rule +r) dQ = w(x), 
Q 2 


OY 4m 4m or 
whence . 
BE es rlth ay (S52 a las 
Thus we have 
(2.9) u(p) =m (h)+%4(p), pin &, 
where ; 
(2.10) %,(p) = ae ce (eee ee (So )has. 


Clearly, v, satisfies the Helmholtz equation in &, since u and u, do so, and re- 
presents the wave function due to layers of simple and double sources on © of 
densities Ou/On and —u, respectively. In exactly analogous fashion we can 
show that 


u(p) = U(h) +%(p), p in &, 


4 Ou [ e—tkR r,) e—tkR 
Uap) = af Ven ( R (on ( ee )}as. 

SOMMERFELD’S radiation condition is frequently used in diffraction problems 
to show that for certain Helmholtz wave functions uw, (f) = 0, so that u(p) =v, (p). 
It is shown below, in § 4, that this is actually equivalent to SOMMERFELD’S con- 
dition for solutions of the Helmholtz equation. Meanwhile, let us use the ter- 
minology that a Helmholtz wave function «(p) (defined in an exterior domain) 
is a radiation function (respectively, antivadiation function) if u,(p) =0 (respecti- 
vely, u,(p) = 0). 

Next, we should like to show that v,(f) is always a radiation function. In 
view of (2.7) this is equivalent to 


where 


(2.11) (s — 8k) rile, r) =0. 


Now, v, =“ —u, whence 3, =% —%, and we know that 


Let us write 


make 5 baht 

ty (&,7) = a 2 ( Vee i 2 (a) 
__ 4 (W, (©) —W, (#)) 4 (e-#7—A) eg Le) lca 
2kyr 2 ky DN) 2 ky aie) 


Then, since u, is an entire function, w,(a,7) is defined for ally and tends to the 
limit u, (a) for 7+0. Consequently, we must have 


Wy (a) = wy(x%) and u(x) =9w, (x) + 3 2(x), 


whence 
W, (#) = W, (x) = u, (x) 
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and : 
(5. — ik) rth, (00, 7) = e7**? w, (aw) =e” my (x). 
Y 

Combining these results we see that (2.11) is indeed satisfied. Let us also notice 
the following lemma. 

Lemma: An entire radiation function u(p) must vanish identically. 

To prove this observe that for an aribtrary p we have u(p) =v,(p) where 
v, (p) is the surface integral (2.10) taken over any closed surface S not containing #, 
since u is entire. On making & shrink to a point we get u(p) =O. 

The lemma shows that any decomposition of a Helmholtz wave function into 
a sum of a radiation function and an entire wave function is unique. Hence, we 
can summarize our results on Helmholtz wave functions as follows. 

Theorem 3: Every Helmholtz wave function u(p) for an exterior domain & has 
a unique decomposition 


u(p) = u,(p) + 4 (2) 


where u,(p) is an entire Helmholtz wave function and v,(p) is a radiation function 
for &. The functions u, and v, are given by (2.7) and (2.10) respectively. 


§ 3. A Radiation Condition for the Diffraction Problem 
The principal result of this section is 


Theorem 4: The diffraction wave function u(p), defined by 


lim e'*! p(p, t) = (6), 


t— oo 


is necessarily a vadiation function. 


For the proof we consider the pulse wave function (Pf, t). If 7)(p) is the 
least radius such that G(p,79(#)) contains the boundary & of ¥ and the set of 


points where /(p, k) = — °PolP: 9+) +0 then §4(p,t,7) is defined for r > 79 (6), 
t>0O and 
(3-1) Gold, t,7) =P EC (rt) 


there. Notice that K(f, t) is defined for all r while G(, 7) is defined for t>7)(). 
Moreover, we know that @o(p,t)=0 outside S(p, 79 (pf) +t), by Theorem 14, 
whence %9(p, 1,7”) =0 for r>7)(p) +4. Thus 


G(p,r +t) =—F(p,r—t) for r>7)(p) +¢, 


or 
G(p,t + 2t) =— F(p,t) for t>7,(p) and t>0. 


This shows that G(p, t) is independent of t for all t>79(f), i.e., on its whole 
domain of definition. Hence, since a constant may be regarded as a function 
of y—t, we can combine G with F in (3.1) and write 


(3.2) Goldy 8) ee ee or 


if 
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Also, since @ (pf, t, 7) =0 for r—t>7,(p), we see that 


(3.3) Fp, t)-= 0.'for, t>16()s 


Let us see that these results imply for the diffraction wave function (p). 
We have 


On integrating this identity over (pf, v) and reversing the order of the integra- 
tions we get 


7) = fet COM Dat, 1) ae = penser) at 
0 : 0 
me a fe FE, ‘dc = = fe CrP (ph, thd 


where at the last step we assume that 7>7,)(p) and use (3.3). On comparing this 
with the general form (2.6) for #(f,7) we see that u,(p)=0, which proves 
Theorem 4. 

Our result suggests that u,(p) =0 is a “‘natural” radiation condition for the 
diffraction problem as formulated in Sect.1. In the next section we compare 
this with other radiation conditions, including SOMMERFELD’s. In particular, 
we shall prove that the conditions 


U, = 0 


characterize the diffraction wave function uniquely, by showing that u,(f) =0 
implies the usual Sommerfeld condition (Theorem 5). 


§ 4. Other Radiation Conditions 
We shall refer to the equation u,(p) =0 as R. C. 1 to distinguish it from the 
other radiation conditions we are about to introduce. A more explicit statement 
of R. C.1 is derivable from (2.7). It is 


R.C.4 [ (252 — iaral as —0 for all p and r>”%(f). 
S(b,7) 
Now, notice that solutions of the Helmholtz equation in & satisfy the identity 


wn f {2 -ttrah as =F [ (Go—sau)as + 
vy? Or (a Y 
S(p,7) S (2,7) 
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Hence, since the last term on the right tends to zero with 1/7 for each p, R.C.4 
implies 
Ree lim 1 [ (Gp-tealas =0 for all p. 
S (6,7) : 
Conversely, since the term on the left equals 4 e—**’u,(p), we see that RieZ 
implies R.C.1. Thus R.C.1 and R.C.2 are equivalent. 
Next, if we apply ScHwarRz’s inequality to the integral in KC. 2,ewe get 


ts [ (Gi —iku)asf sam f |Z —ikulas. 
Y or ae ee oY 


(P,1) 

Thus R.C.2 is implied by 

R.C.3 lim [ |E-ttafias =0 for all p. 
eaten 

Since 


Ou(a+rr) 2 lhe 
oe iku(w +r)}| aQ 


i a —ikul dS = f 
OY 

GS (a,7) Q 
we see that R.C.3 is implied by 
RoC lim (3 — ik) =0 for all x, uniformly in r. 

%—>0O Y 

This is SOMMERFELD’S radiation condition. Finally, we state a weakened version 
of R.C.3: 


Rees lm if ea ~ ikul aS =0 for a single center O. 
%—> 0O or 
S(O,7) 

The author has found no discussion in the literature about whether SOMMER- 
FELD’S condition is to hold for only one, or more than one, center. However, 
in deriving results like the representation theorem, representing w in terms of 
a surface integral (2.10), it is always tacitly assumed to hold for all centers 
[1, 9, 14]. On the other hand, RELLICH proved his uniqueness theorem [7] using 
only R.C.5, which suggests that R.C.5 may imply the stronger conditions. We 
verify this by proving. 


Theorem 5: Let u(p) satisfy the Helmholtz equation in an exterioy domain &. 
Then Radiation Conditions 1 through 5 are equivalent for uw; i.e., each implies all 
the others. 


Proof: We have seen that R.C.1 and R.C. 2 are equivalent and that R.C. 3 
and R.C. 4 imply R.C. 2 and R.C.1. Moreover, R.C.4 implies the representation 
theorem and hence the expansion theorem [/4] which states that for spherical 
coordinates (7, 7) relative to an arbitrary origin, 


etkr 


u(p) = — ee ine rece 
n=0 


ia 


where the series converges uniformly and can be differentiated term by term 
any number of times, the resulting series all converging uniformly. This makes 
it clear that R.C.1 implies conditions R.C. 2 through R.C.5. It remains to 
prove that R.C. 5 implies R.C.1. 
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We must show that R.C.5 implies #(f)=0. To do this let us use (2.8), 
taking for SG’ the sphere G(O,7) of R.C.5. Then n=r and so 


a ikR a / eikR 
24 (2) = J ter 4 Bierce )}as 
S(O 
ikR : ikR , 
= [te g —ulik #) = — las = [7r(& iku) x 
2 


Y ikR . eles, gikR aA OR gikR 
X ee ee ae aE? goat sag dQ, 
2 


where the last integrals are with respect to solid angle 2 viewed from O. Let us 
write 


4nu,(p) = +1,+ 1s, 


where J,, J,, and J; are the last three integrals, and estimate these terms separa- 
tely. Notice that 
R2=7—2 or <1+ 4 


R 
where d is the distance from # to O. Thus - <2 for R2d. On using this estimate 
and SCHWARz’s inequality we get 


HS Ai cemtd! aS, 


whence lim J,=0, by R.C. 5. By the same reasoning 
1 © 


2 
|Z,|2< se | Weliase — 
2 


Now R.C.5 implies {\uPae= O(% =) [14]. Moreover, lim a =1 uniformly 


on Q, as is easily verified. Thus lim J,=0. Finally, since 


|Z]? S 64x [ |ul2a2 —O(7,), by R.C.5 [14], 


whence lim J;=0. Collecting these estimates we have lim ~,(p) =0. Thus, since 
%—>0o 1%—> CO 


u,(p) is independent of 7, “,(p)=0, which completes the proof. 


§ 5. The Vector Diffraction Problem 


The initial-boundary value problem for vector, or electromagnetic, waves 
with time-harmonic sources is to find the vector fields E(p, t) and H(f, t) uniquely 
characterized by the conditions 


Vxbt = =0, VxXH- 2 —e NG), pinB,1>0) 


E(f,0+)=0, H(p,0+)=0, pin&, 


RCE G1) =O" son Spt 2:0" 
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Here n is a normal vector to G, and we assume that F(p, k) vanishes outside a 
sphere. The corresponding pulse solution is defined by 


VxE,+—2=0, VxH——~=-0, pin, t>0, 


E,(?,0+) =—F(¢,k), Ho(o,0+)=0, pins, 
nxE,(~,t) =0 on SG, 420. 
The problems are related by 
t 


E(p,t) =f E,(p,t—t)e***dt, 


0 
H(f, t) =| Hy(p.t— Te ans 
as in the scalar case, and if the integrals 
fet E, (p,t)dt and fe™ (p, t) dt 
exist then : ‘ 
E(p,t) =A(p,t) +e" A(p), 
H(p,t) = B(p,t) +e Bip) 


where A(, t) and B(pf, t) are transient and 
A(p) = fe Ey(p, t) dt, 
0 


B(p) = fe Hy(p, 0) dt. 
0 
Finally, we can verify that 
VxAttRB=F, VxB—7tkA=0, Yin 3, 


C4) nx A(p)=0 on G, 


but that these conditions fail to determine the fields uniquely. 


A radiation condition that provides a unique solution to (5.1) was discovered 
by S. SILveR [8] and, independently, by C. MULLER [6]. It is 


lim {rx (VxA) +7kA}=0, uniformly in fr. 


We shall show that the diffraction field A defined above must satisfy the Silver- 
Miiller condition. 


If S(p,79(p)) contains G and the set on which F(f, k) +0 then u=V-E, 
and v=V- Hy satisfy du/0t = 0v/dt =0 and u(p, 0+) =v(p, 0+) =0 at points 
outside &(p,79(f)). Thus V- E,(p,t) =V- Hy(p, ) =0 outside S (p, 7%9(p)) for 
t>0O and therefore 


2E 
AK, — a = 


E,(6,0-+) =0, SE OF) 9 


outside S(p, 79(f)). 
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It follows by Theorem 1 that E)(p, t)=0 outside G(f, 79(p) +2). This implies 


R.C.4 [ oP -ieralas =o for all p and r >7,(), 
S(?,7) 


exactly as in the scalar case. To relate this to the Silver-Miiller condition we 
use the identity 


rx (VxA) =V(r- A) — oO r=Trr, 


valid for all fields A of Class C1, which is easily verified by using rectangular 
coordinates. This gives 


r {ex (PX A) +ikAY=V(r- A) — {2¢4) = ikr Al. 
Notice that on identifying # with its position vector a, and the integration point 
g with #)+ 7, so that A(q) =A(r+a,), we have 
Vir-A) =A+D/,(r-A), 
where Vy refers to 2. Moreover 
f Yor A) dQ = arn je (Vx B)as| =O 
2 


S (6,7) 
by MAXwWELL’s equations (5.1) and STOKEs’ theorem. Thus 


x Or 
S (2,7) S (2,7) 


1 (ee —ikralas=—4 [ @xoxa) + ik A} dS + A(p,7) 


where 


It follows that R.C.1 is equivalent to 


R.C.2 lim 1 a ex (Vx A) +ikA}dS =0 for all f. 
aes Cc (2,7) 

The remainder of the discussion parallels that for the scalar case and the details 

are omitted. As a final result we can state the following analogue of Theorem 5. 


(Cf. also [15)). 
Theorem 6: Let A(p) satisfy the vector Helmholtz equation 


Vx (VxA)—PA=0 


in an exterior domain &. Then conditions R.C.1, R.C. 2, the Silver-Miiller condi- 
tion and 


R.C.3 lim f |rx(VxA) +ikAl?dS =0 


%¥—> 0O S(0,7) 


for a single center O, are all equivalent for A. 
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Einleitung 


In den Jahren 1908 bis 1920 hat P. KOEBE in einer Reihe von Arbeiten das 
Problem behandelt, ein m-fach zusammenhangendes Gebiet G mit oo€G auf ein 
nur von Kreisen berandetes, sog. Vollkreisgebiet V konform und in oo normiert 
abzubilden. Zur Losung dieser Aufgabe schlagt er zwei Iterationsverfahren vor, 
das ,,iterierende Verfahren“ und das ,,Iterationsverfahren“, deren Hauptgedanke 
es ist, die gesuchte Abbildung durch Hintereinanderschaltung einer Folge von 
konformen Abbildungen einfach zusammenhangender Gebiete zu gewinnen. Wir 
wollen diese Verfahren funktionentheoretische Iterationsverfahren nennen, im 
Unterschied zu Methoden, die z.B. im Fall 7 =2 auf die Lésung eines Paares 
von Integralgleichungen abzielen ([7]; [8], S. 137—147). 

Da es KoEBE in erster Linie auf die Existenz der Vollkreisabbildung ankam, 
konnte er sich mit einem bloBen Konvergenzbeweis fiir seine Iterationsverfahren 
begniigen. Die Geschwindigkeit der Konvergenz und die Frage der gleichmaBigen 
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Konvergenz in ganz G blieben demgemaéB auBer Betracht. Vom praktischen 
Standpunkt aus sind jedoch solche Fehlerabschatzungen erwiinscht, weil sie ein 
Urteil tiber die Giite des Verfahrens erlauben, selbst dann, wenn jeder Iterations- 
schritt wieder nur angendhert ausgefiihrt werden kann, wie das in der Praxis 
der Fall ist}. 

In Kapitel III dieser Arbeit soll nun gezeigt werden, wie die Grundidee 
Koesks, die sich durch das Stichwort ,,Erhéhung der Spiegelungsfahigkeit der 
Bilder von G‘‘ markieren laBt, beweistechnisch so verfeinert werden kann, daB 
man zu brauchbaren Fehlerabschatzungen kommt. Wichtige Hilfsmittel dazu 
sind die Satze 1.1 und 1.2 von Kapitel I. AuB8erdem sollen in Kapitel III die zu 
den Koebeschen analogen Verfahren bei Normierung im endlichen auf ihre Kon- 
vergenzgeschwindigkeit untersucht werden (§ 9, B und C), wobei sich im einen 
Falle eine Verbesserung fritherer Abschatzungen [5] zum Verfahren von KOMATU 
ergibt. Fiir einige der in Kapitel III angegebenen Iterationsverfahren sind Arbei- 
ten im Gange mit dem Ziel, an einfachen Beispielen die Giite der gewonnenen 
Abschatzungen nachzupriifen; vgl. dazu auch § 10. 

Unabhangig davon betrachten wir in Kapitel II hauptsachlich konforme Ab- 
bildungen nahezu kreisringférmiger Gebiete auf Kreisringe. Darin wird die Aus- 
sage prazisiert: Ist G ein nahezu kreisringférmiges Gebiet und f(z) die (normierte) 
konforme Abbildung dieses Kreisrings auf G, so ist /(z) ,,nahe bei‘‘ der Identitat z. 
Ein erstes Ergebnis zu diesem Thema [8] wird dabei verbessert. In § 5 findet 
sich weiter eine Aussage tiber das Verhalten der (normierten) Abbildungsfunktion 
eines Gebiets, das von |w|<1 ,,wenig abweicht‘*. Die Abweichung findet dabei 
nicht wie gewohnlich bei |w|=1 statt; G wird sogar von |w|=1 berandet, aber 
auBerdem noch von einem in |w|<@ verlaufenden Randkontinuum (g klein), 
so daB G zwar ein Ringgebiet ist, aber doch von |w|<1 ,,wenig abweicht“. 


I. Vorbereitende Satze 


Als Vorbereitung fiir unsere spaiteren Untersuchungen wollen wir zundchst 
einige Satze zusammenstellen, deren Inhalt zum gréBten Teil bekannt ist. 


§ 1. Hilfssdtze fiir Ringgebiete 


A. Wir beginnen mit der Behandlung zweier Extremalprobleme fiir endliche 
Ringgebiete. 


Satz1.1. Es set G ein endliches Ringgebiet vom Modul M(>1) mit der inneren 
bzw. duperen Randkomponente C, bzw. C,. Bezeichnet dann F (C,;) den von C; um- 
schlossenen Fldcheninhalt (i =1, 2), so gilt 


Ly 1 Ly 
(1.1) Fehr 


A terber tritt das Gleichheitszeichen genau dann ein, wenn G ein konzentrischer Kreis- 
ring mit dem Radienverhdlinis M ist. 


' ; Eine Vollkreisabbildung eines zweifach zusammenhangenden Gebiets nach dem 
iterierenden Verfahren hat Léscu in [19] durchgefiihrt. 

2 Falls der Flacheninhalt nicht existiert, so ist fiir 7—1 der duBere, fiir i—2 der 
innere Peano-Jordan-Inhalt zu nehmen. 


Durchfiihrung der konformen Abbildung 154 


Der Satz stammt von CARLEMAN ([3], S. 212) und wird am besten tiber einen 
Flachensatz bewiesen; vgl. dazu [10], S.177 und [9]. In [27], S. 192 wird (1.1) 
mit Hilfe des Prinzips der zirkularen Symmetrisierung gewonnen, und schlieB- 
lich sei noch bemerkt, daB Satz 1.1 einen Teil der Ergebnisse in [4] enthalt, da 
dort C, als Kreis angenommen ist. 

Etwas tiefer liegend ist folgender 


Satz 1.2. Es bezeichne D(C,) den Durchmesser der Randkomponente C, von G 
aus Satz 1.1. Dann gilt 
(12) [DG)PSKM)-F(C) mit KM) =" Y Tea 

n ungerade 

Uber die Annahme des Gleichheitszeichens laBt sich sagen: Wird 1< |z| <M? 
durch w=f,(z) auf einen zur u-Achse (w=u-+iv) symmetrischen und ldngs 
—1SuS +1 geschlitzten Parallelstreifen konform abgebildet, so tritt in (1.2) sicher 
dann Gleichheit ein, wenn G von Cy: —1Su<+1 sowie von Cy: fy(\z| =M) be- 
vandet wird. Jedes weitere Extremalgebiet ergibt sich daraus durch eine ganze 
lineare Abbildung. 

Diesen Satz 1.2 hatten wir in einer vorbereitenden Arbeit [6] mit Hilfe des 
Prinzips der zirkularen Symmetrisierung und Saétzen von AHLFORS und SCHIFFER 


bewiesen. Dabei laBt sich K(J/Z) geschlossen angeben: K(M) = —= log k’, wenn 
k’ den zu k =k(M~?) komplementaren Modul bezeichnet. Das asymptotische Ver- 
halten von K(M) laBt sich durch 

BAM oe eo E03) 
und S 1 , 

K(M) = Se + O(M --1) (M—>1-+) 
charakterisieren; [6] enthalt auBerdem eine kleine Tabelle fiir K(M). Schlieb- 
lich bringt eine direkte Abschatzung von [D(C,)]?/F(C,) mit gelaufigen length- 
area-Methoden K(M)< oF (M>1); die Abschatzung ist fiir nahe bei 1 


2 - log 
liegende Werte von JM recht scharf. 


Fiir unsere Anwendungen von Satz 1.2 im Kapitel III merken wir noch an: 
Fiir wachsendes M fallt K(M) monoton. Dies ergibt sich aus der Reihendarstellung 
fiir K(M) sofort. 


B. Mit dem Fragenkreis des ,,nachsten Randpunkts“ bei Ringgebieten hangt 
eng zusammen 


Satz 1.3. Es set G ein endliches Ringgebiet vom Modul -. (>1), welches auBen 
von |w|=1 berandet wird und das w =O nicht enthalt. Dann legen alle Punkte 


der inneren Randkomponente C, von G im Kreis |ol< a falls rS> ist. 


Beweis. Bekanntlich? wird der Abstand von Punkten von C, vom Null- 
punkt dann maximal, wenn C, eine Strecke 0S w<A(r) ist, so daB A(7) abzu- 
schatzen bleibt. Nun kann man aber das von |w| =1 und dieser Strecke berandete 


3 Vel. dazu das ,,Problem des nachsten Randpunkts‘“ in [7/7] sowie bei TEICH- 
MULLER [23], insbesondere S. 631—633. 
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Gebiet dadurch erhalten, daB man das von |w|=1 und —L(7)SoS+L(7) 


berandete Gebiet der linearen Abbildung w= ot unterwirft, wobei dann 
ee LAG . 
A(r)=4573 ) wird. Wegen 
7 1+78" \2 
LO) 27, I (y=) 2r 
20], S.293—295) ist daher Aye da g(x) = 2" in, O55 % 524 
(20), 8.293 2 Ae < ft, in 054s 


monoton wachst. — Da r<|w|<1 den Modul 1/7 hat, sieht man, daB stets 


A(r)>r ist; A(r)<—*" 


PaNTE gilt also fiir nahe bei 1 liegende 7 nicht mehr. 
if 


§ 2. Eine Abschatzung von SPECHT 


In §4 werden wir auf folgende spezielle Abschatzung zuriickgreifen, die 
SPECHT ([22], S.187—188) im Verlaufe der Untersuchung konformer Abbildungen 
nahezu kreisférmiger Gebiete bewiesen hat. 


Satz 2.14. Die Funktion o0(g) sei mit 2 periodisch, in <0, 27> absolut stetig, 
und es gelte 


o(y)20>0 sowie 


o'(@) aK ¥ lle o. 
cols fir fast alle p 


Die Funktion y(t) sei in (0, 2a) absolut stetig und monoton wachsend, und p(t) —t 
set mit 2% periodisch. Dann gilt fiir fast alle } 


2% 
1 O—1 
(2.4) Li, foselo) ctg —— dt 
0 


2a 
<K —-/ by weak, 
SK |2l0g2 rss (p(t) —#) ctg 5 atl. 
0 


Der Beweis erfolgt durch zweimalige partielle Integration. 


II. Uber die konforme Abbildung nahezu kreisringférmiger und nahezu 
kreisformiger Gebiete 

Wird ein mit |w|<14 nahezu tibereinstimmendes Gebiet G durch w =/(z) kon- 
form und normiert auf |z|<1 abgebildet, so ist zu erwarten, daB (unter genaueren 
Annahmen) | f(z) —z| klein ist. Einen entsprechenden Satz fiir Ringgebiete hatten 
wir schon in [5], S. 873—877 bewiesen, doch lieB die dort durchgefiihrte Reduk- 
tion auf den Fall einfachen Zusammenhangs nur die GréBenordnung von [#(z) —2| 
ohne brauchbare Konstanten abschatzen. Das Problem soll deshalb erneut auf- 
gegriffen werden, obwohl dies fiir Kapitel III nicht erforderlich ware, und es 
wird sich zeigen, da8 wir unter Bentitzung von Zusammenhingen konjugierter 


Funktionen im Kreisring auf direktem Wege zu brauchbaren Abschatzungen 
fiir | f(z) —z| kommen kénnen. 


§ 3. Zusammenhang konjugierter Funktionen im Kreisring 
Gegeben sei eine in R<|z|<1 regulare und in RX|z|<1 stetige Funktion 


i (2) =u(re’®) +iv(re). Wir benotigen zunachst verschiedene Zusammenhange 
zwischen den Funktionen u und v auf dem Rand unseres Kreisrings. 
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Satz 3.1. Ist f(z) in R<|z|<1 regulary, in RS|z|<1 stetig, so gilt 


20 
(3.4) He") = 40) +5 fie)cte2= dt (0S S22), 
wober 
1 40 d R 
3.2) A(d) = gai JO pas “und also | A(8)|< = max | (2) 
gilt. 


Eine Zerlegung von (3.1) in Real- und Ima- 


ginarteil ergibt 
20 


(ta) u(e®) =mA(O) +S" fo(e cg Eas, 


0 
2 


(3.1b) v(ei*) = 346) +50 [ ule’ gets 


0 


samtliche auch im folgenden auftretenden un- 
eigentlichen Integrale sind als Cauchysche Haupt- 
werte zu nehmen. 

Entsprechend lautet der Zusammenhang zwischen w und v auf |z| =R 


Satz 3.1’. Ist f(z) in R<|z|<1 regular, in R<|z|<1 stetig, so gilt 


Fig. 1 


64) 9 7(Re*%) = BOO) + ae f Hel cg 2h dt = S0 S20), 


wober 


, at gece ONES dz ik 
(3.2') Bd) = oa fio ON eta und also | B(#)|< Tejas | #(z)| 
galt. 
Eine Zerlegung von (3.1) in Real- und Imaginarteil ergibt 
22 
(3.1’a) u(Re'*) =RB(d) +a f (Re) eed 
0 
22 
(3.1'b) v(Re!*) = 3 B(d) +5 fuRe's otg Oat. 


tt) 
Man erkennt, daB z.B. durch (3.1a) die zu v 1m Ring konjugierte Funktion u 
in Zusammenhang gebracht wird mit der zu v in |z|<1 konjugierten Funktion, 
mit einer durch #A(#) ausgedriickten Korrektur. 
Wir beweisen etwa die Beziehung (3.1). Mit den Bezeichnungen von Fig. 1 
gilt ys fie ete Azam: ni {16 eet dz 
Dt 


oa ere_z gz ee_z z 


> f (ei) 4 +g fine eit yotgtea di (o> 0). 


Daraus folgt (3.1) unmittelbar. Der Beweis von (3.1’) verlauft entsprechend. 
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Die eben abgeleiteten Formeln (3.1) und (3.1’) haben noch den Nachteil, daB 
die in ihnen vorkommenden Korrekturen A(#), B(#) noch von u und v abhangen. 
Wir leiten nun zwei zu (3.1b) und (3.1°b) entsprechende Beziehungen ab, in 
welchen sich die Korrektur allein durch wu ausdriickt. 

Satz 3.2. Ist f(z) in R<|z|<1 regular, in RS|2|S1 stetig, so gilt mit einer 
Konstanten Bo 


steps Sipe ne het) — rete sont 
(3.3) ce 
oan ctg a, 


0 


co 
e*) = B+ >, woes . fem —u(Re'')\.sinn(t — 8) dt + 
(3.4) rie 2a 
* aq f mle cg de 
0 
Zum Beweis kniipfen wir an eine zunachst fiir R<7<1 giiltige Darstellungs- 
formel von ViILLaT [24] an. Setzt man 


27 22 
a, = [ we") cos nt dt, a= + [ u(Re') cos nt dt 


2a 27 eet): 


b= = f ule") sin nt dt, b= — f u(Rel sin nt dt 
0 0 


so gilt fir R<r<1 


~ (olite fee 
v( FD " = ates cosn ? — a, sinn #) — 
~ 2m aoe (b;, cosa 9 — aj, sinn #). 
Wir untersuchen diese Relation bei festem # fiir y+1—. Eine Anwendung des 


Abelschen Stetigkeitssatzes zeigt, daB die zweite Reihe den Grenzwert 


co 
4 ‘ 
2.) —, (b, cos 2 9 — aj, sin n 9) 


besitzt. Um die erste Reihe zu beurteilen, schreiben wir 


LR)E SCT ig to ee , (Riv)" + (7 RY” 
SRS Rt~ R-” , 


so daB der vom zweiten Term herriihrende Beitrag den Grenzwert 


Tae’ : 
ace | cosn¥ — a, sinn ®) 
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hat. SchlieBlich bleibt noch >’ (,, cosn#—a, sin n 0) ”” fiir y—>1— zu beurteilen, 


n=1 
also der Abel-Limes der zu >’ (a, cosnt-+},,sinnt) konjugierten Fourier-Reihe. 


n=1 


Da auch v(z) in R<|z|<1 stetig ist, existiert derselbe, und sein Wert ist 


Fassen wir zusammen, so ergibt sich 


v(e*) = 844 y ni (6, cosnd — a, sinn}) — 


n=1 
27 


[o.e) 
2 , aa ; —t 
roe (b,,cosn9 — a, sinn®) tae fue ctg = ai, 


0 


also gerade (3.4). Der Beweis von (3.3) verlauft entsprechend. 


$4. Konforme Abbildung nahezu kreisringformiger Gebrete 
Nunmehr sind wir in der Lage, unseren angekiindigten Satz zu beweisen. 


Satz 4.1. Das Gebiet G der w-Ebene sei von zwei beziiglich w=0 sternigen 
Jordan-Kurven C, und C, berandet; thre Polarkoordinatendarstellungen seven 0, (p) 


und 02(9). 
a) Die Funktionen 0;(p) seen in (0, 27> absolut stetig. 


b) Fur ein e>0 gelte 


<a(p)<1+s und (b) ~t- <i <ite (0<pSs2za). 


all pee 


c) Fiir die nach a) fast iiberall existierenden Ableitungen 9; (p) gelte 


(Cy) und — (Cs) 


das heiBt es ist 


ltgu(9)|<N-2e baw. |tgx,(y)|<N-e, 


wenn x;(p) den fast tiberall existierenden Winkel zwischen |w| =0,(y) und C; im 
Punkte (p; 0;(p)) bezeichnet. 

Bildet dann w=F(z) den Ring R<|z|<1 konform auf G ab, wober F(R) >0 
ist, so gibt es eine nur von R und N abhingige Konstante K, so da gult 


(4.1) |F(z) —2|SK-e (Rie 2) 54). 


Der Beweis wird zeigen, daB man im Falle «<1 hierbei 


(42)  K=[1+2N(2log2 +1)] +7 25a5 + 2N- 4142 4 4 20) 
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wahlen darf; ist speziell 0.(y)=F fiir alle y, so gentgt 
(4.2') =[1 1+ N(2log2 +4)] + Sat (+N), 


so daB diese Konstante fiir R->0 in die von Specut fiir den Fall einfachen 
Zusammenhangs angegebene Zah] 1+ N(2log 2 +4) tibergeht. 


Zum Beweis unseres Satzes betrachten wir die Funktion 
f(e) =log2@) = uz) +iv@)  (RS|z|S1); 


F(R) 
R 


durch die Forderung v (Rk) = arg ( 
Dabei ist wegen b) 


| ( 2)| = |log |=! 


)=0 ist die Funktion eindeutig festgelegt. 


|| s tog (i 2) te ~auf|2| 4) und .|z) = 


und unser nichstes Ziel ist es, auch v(z) auf |z|=4 und |z|=R abzuschatzen. 


a) Abschétzung von |v(Re'*)|. Verwenden wir (3.3) und schatzen zuerst die 
Reihe ab, so ergibt sich 


Der letzte Term in (3.3) gibt AnlaB zu 


(4.3) 


2a 
1 O—t 
aq [8 e2((d)) ete 7 a 
0 


22 
< np ot Gata 
 Ne|2I0g2 aa | (7 t) ctg 5 atl ; 
0 


dabei wurde Satz 2.1 verwendet, dessen Voraussetzungen hier alle erfiillt sind. 
Zur Abschatzung der rechten Seite von (4.3) verwenden wir nun (3.1’a) und 
erhalten mit (3.2’) 

2% 
Sait v (Re) ctg ae 
27 


0 


=a a (s + 22); 


dabei haben wir |v (e**)| grob durch 2% abgeschatzt. Insgesamt ergibt sich so 


lo(Re*) —6,| S4e_—"*_ Ne 


U—RpP “ 2log 2 +2 + +t% (e+ 2m)|=a-e. 


Wegen v(R)=0 entnimmt man daraus |f)|<«-«, also (fiir fast alle # und 
wegen der Stetigkeit von v(Re*®) sogar fiir alle d) 


(4.4) |v(Re)| < 2a-e. 


b) Abschéitzung von |v (e**)|. Aus (3.4) folgt wie oben sofort 
22 
1 


i O—t 
= f «lc ‘) ctg dt : 


0 


|v(e*)|Sa-e+4e- 
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Hierin ist der letzte Term nach Satz 2.1 und (3.1a) 


2n 
= salstey 2 Ong 
= e|210g2 sa | (7 ) ete oa <Ne 


0. 


Es ergibt sich somit (wieder fiir alle 8€ <0, 22>) 


|v (2) | <a- et {sera 


also ist 


+Ne e|2log2 ++ ae (e+2n)]| =2 


etsy: 
iz) arto es fire lz und! slz| = 4 . 


Daraus folgt nun schlieBlich 


|F(e) —2| <|72 — 1] <||7@ arg AE se 4 20, 


~4|4 


o Boe 1+R 
(1—R)? NTR 


|F) —2| Se|1 + 2N(e + 210g 2) 4 (e+ 22)|, 


giiltig fiir |z| =1 und |z| =R, also in RS|z| <1. 
Fihrt man den Beweis fiir den Sonderfall 9,(~)=R durch, der in [5] be- 
handelt ist, so erhalt man zunachst aus (3.3) entsprechend zu oben 


/D pf i 8R 
{Re =|0(Re*| Sees 
und (3.1b) liefert weiter 
2x 
ee 1+R q 8R ‘ 1 th it O—t 
lu(é*)|s =e GER! Seale u(e ) ctg —— at}, 


0 
wovon der letzte Term nach Satz 2.1 abgeschatzt wird und 


2% 
1 it o—t 
rene 
Ne |210g2— f ole) ctg dt 


ist. Mit (3.1a) ergibt sich 


jo(e9)| Se | PA, -8R +N. 2log2+Ne+ oto, 8RNe| 


S e|N(e + 2log 2) += Saeed +Ne], 


(1—R)8 


also 


|F(2) —2| Se/1+N(e + 210g 2) ! ae (44 Ne). 


Daraus folgt die Aussage (4.2’) fiir die Konstante K. 


§5. Konforme Abbildung nahezu kreisformiger Gebrete 


Zum Beweis der Satze 9.3 und 9.4 bendtigen wir noch einen Satz, der sich 
mit der konformen Abbildung von Ringgebieten beschaftigt, bei denen eine Rand- 
komponente klein ist, so daB die Abbildungsfunktion dort , nahezu“ eine heb- 
bare Singularitat besitzt. 
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Satz 5.1. Gegeben sei ein zweifach zusammenhdngendes, endliches Gebtet G 
der w-Ebene, welches auBen von | «| —4 berandet wird und das w =O nicht enthalt. 
Es bezeichne w=f(z) die durch die Bedingung {(1)=4 normierte konforme Ab- 
bildung von r<|z|<1 auf G. Dann gilt fir eine absolute Konstante C 


[f(2) —2|SC-r (7 Siz) 1 O71): 


die kleinste derartige Konstante Cy gentigt 4C)S12,6. 


Beweis. Wir zerlegen zunachst f(z) in eine in |z|< - regulare Funktion F(z) 


; : i har 
und eine Korrekturfunktion k(z). Fiir jedes z mit r<|z|< a gilt 


AGT Roeee f(o)/6 ee ay ek ES ( +) 
Z. 220i ee os 201 f—z te 4<<|2]|<n< ; 
[ln [el=re 
also fiir f(z) die Zerlegung f(z) =k(z) + F(z), wobei F(z) in |z|< “ regular ist 
und F(0) =0 ist. 
Fir die in |z|>7 regulare Korrektur k(z) ist nun eine méglichst gute Ab- 
schatzung wichtig. Eine erste grébere Abschatzung ergibt sich sofort mit 
Satz 1.3; fiir jedes 7, mit r<7,<1 gilt namlich 


4 1 
Oe ar eer (l2| =), 


und 7,7 liefert 


(5-1) AQ st =a) (lz|=1). 


Wir merken gleich an, da8 ftir O<7<+ die Funktionen se und c, (7) monoton 


zunehmen, und daf c,(%) <1 ist. Fiir eine etwas genauere Beurteilung von |k (2) | 
schreiben wir 


(65.2) bey = f1@|-4,-|@- = ipiar dt—*. ge dt, 


[lan [S|=ry [Sl=n 


wo C eine beliebige Konstante sein darf. Beriicksichtigen wir, daB die Menge 
der Punkte /(¢) (|¢|=7,) einen Durchmesser <2L(7,) <4”, hat‘, so gibt es nach 
einem elementargeometrischen Satz von JuNG® einen Kreis vom Radius 4n,/V3, 


in welchem diese Punktmenge liegt. Nehmen wir C als dessen Mittelpunkt, so 
ergibt sich aus (5.2) 


47 4% 


| k(z) | = ae Tan (ae 
und 7,  liefert 
472 47 
(5-3) |k(z)|S Bizi—n | 1+472 (J2|>7). 


* Dies geht aus den Ergebnissen von GrorzscH in [12] hervor; L(v) ist die im 
Beweis von Satz 1.3 auftretende GréRe. 
SWAN, 9,18, IST, Se TEXO: 
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Es sei nun weiter 7 auf O<7<3 beschrankt, so daB |k(z)|<1 wird fiir |z|=1. 
Da aber | f(z)|>1 ist (|z|=1), haben wir A arg/(z)|,,,— ea ae Aarg F(z 


und F(z) hat daher z=0 als einzige Nullstelle in |z|< -. Also ist 


Vlete =1/r> 


H(z) = log ze) =u(z)+7v(z) in |a| <= 
regular erklarbar. Wir merken gleich an: 


(5.4) u(0) =log|F’(0)|S0, da |F’(0)| = 


a4 
ae [4 : hats 


[S]=1 


ist. AuBerdem ist auf |z|=1 unter Beniitzung von (5.1) 


(5.5) | «(2)| = [log |F@)|| <|log(t4 —a”))|=e2(r) — (\2| =1) 
und also fiir jedes 9 mit r<o<1 


ie 1) ule ag] so cate) =e) 


0 


(5.6) | H(z) — H(0)| = 


Diese Abschatzung schreiben wir jetzt um in eine fiir H(z*), wobei * jeweils die 
beziiglich des Einheitskreises vorgenommene Spiegelung bedeuten soll. Es ist 
(also fiir |z| =e, 7<o<1) 


Hi(z®) = log ay) = log {ie ) 1 aH |} = log He + log f Abts) 


wobei der erste Term 


log 2— =1 
psa aoe 


der zweite Term hingegen betraglich 


a (ee AC al aar aa 
an mass a) 


(z)|S 42 und (5.3) verwendet. Setzen 


ist; dabei haben wir [4g 


meee ie 
wir «= —H(0), so folgt aus a 6) 


|H(z*) — «| S|H(2) —H(0)| +A() S 


oder 
| A(z) 


r) + h(o) (lz = -) 


LaBt man hier bei festem z 9—> 7 riicken, so erhalt man 


— 


| (2) — a] 2” g(r) +h(7) =05(7) (lel <<), 
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und dabei ist 


ht) = jog ee ene rs cera 


fiir unsere betrachteten y-Werte monoton zunehmend. Auch c3(7) wachst daher 
monoton mit 7 (0S7S §). 
Wir erhalten nun weiter 
1 
|H(2) — H(0)| $|H@) —a| +|H(0)—a| S2e(7) (lel < +), 
also nach dem Schwarzschen Lemma sogar 


| H(z) — H(0)| S 2¢5(r) r- |2| (|z| <<). 


Diese Ungleichung verwenden wir fiir |z| =4 


(5.7) | H(z) — H(0)| S 27 cg(r) (|2| =4)- 
Es ist namlich unter Beriicksichtigung von (5.4) 
| ©) =e et) = Mages Sale ioe | f [| = (e217?) —¥ 1) (9) cs 027 es(") zn 


F(z) SSFO)| Ss e2realr)__ 4 
Zz == , 


also 


10) Pose 1440) =40) (lel =1). 


2 
Die Normierungsvorschrift f(4) =1 bringt 
|1—F'(0)|<eq(r) 


ein, und somit erhalt man schlieBlich 


\/(@) —2| =| —4) < 2c407) (|2| =1), 
das heiBt 
(5.8) li) —2| 2/4 4), (\z] <1). 


Hierin ist der Koeffizient von 7 fiir unsere r-Werte 5 Bay >7, und da 
if 


trivialerweise oe |/(z) —z|S5r ist, ist (5.8) sogar fiir alle z mit r<|z|S1 
2|>r+ 
giltig. Der Koeffizient von 7 in (5.8) ist auBerdem in Or monoton zu- 


nehmend, und fiir 7 = a3 erhalt man durch Rechnung 


(5.9) | f(z) —z| S12,6-7 (r20)2|s4)0 


5 <v<1 hat man aber 


? 


so daB dies fiir alle ry in 0<7r< _1_ gilt. Fiir die 7 mit 
stets os 


/f(@)—2)S2=4-7<12,6-7 (ry <|z|S1), 


also gilt die Behauptung von Satz 5.4 mit der Konstanten C =12,6 fiir alle 7 
aN OS? il. 
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Um schlieBlich zu sehen, daB C,=4 sein muB, betrachten wir die Abbildung 
von r<|z|<1 durch die lineare Funktion 
— 1+ir z-ir 
1(2) 1—i7- 14+irz° 
Das Bildgebiet erfiillt die Voraussetzungen von Satz 5.1 und es ist /(1) =1. 
Bemerken wir noch 
1f(=1) —(—4)| = 47 +0(77) (7-0), 


so sieht man, daB in unserem Satz keine Konstante C< 4 genommen werden 
kann. Damit ist Satz 5.1 vollstandig bewiesen. — 

AbschlieBend méchten wir noch bemerken, da8B man in Analogie zum Vor- 
gehen bei BIEBERBACH [2], also mit Hilfe von Flachensdtzen, auch zu einer 
Abschatzung von |f(z)—z| kommen kann. Allerdings ergibt sich dann nur 


[f(2) -2| =O(7). 


III. Funktionentheoretische Iterationsverfahren zur konformen Abbildung 
mehrfach zusammenhdangender Gebiete 


Wir kommen nun zum Hauptthema der Arbeit, der Gewinnung von Fehler- 
abschatzungen fiir Iterationsverfahren zur konformen Abbildung mehrfach zu- 
sammenhangender Gebiete. In §7, §8 und §9, A ist G ein n-fach zusammen- 
hangendes Gebiet der z-Ebene mit «€G, welches konform und in oo normiert 
auf ein Vollkreisgebiet V abgebildet werden soll. Eine Abbildung durch /(z) heiBt 
dabei in oo normiert, wenn in der Umgebung von o gilt f(z) =z+ a. el 3 ae 
Hintereinanderschalten normierter Abbildungen gibt wieder eine solche. 


Den Koebeschen Iterationsverfahren liegt folgende Uberlegung zugrunde. Das 
Vollkreisgebiet V laBt sich iiber seine Rander hinaus beliebig oft spiegeln; um- 
gekehrt ist jedes Gebiet, das diese Eigenschaft besitzt (im Sinne von § 6), ein 
Vollkreisgebiet. Daher werden solche Funktionen /,,(z), fiir die /,,(@) an dessen 
Randern oft gespiegelt werden kann, eine gute Approximation fiir die exakte 
Abbildung von G auf V darstellen. Zu den /,,(z) kommt KOEBE tiber das ,,ite- 
rierende Verfahren“ (§ 7) und das ,,Iterationsverfahren™ (§ 8). 

Fehlerabschatzungen waren bisher fiir diese Verfahren nicht bekannt. Auf 
Anregung des Verf. hat neuerdings HUBNER [/4] in dieser Richtung Fortschritte 
erzielt und Aussagen entsprechend den Satzen 9.1 und 9.2, aber mit schwacherer 
Konvergenz, gewonnen; die Beweisfiihrung bei HUBNER ist allerdings auf den 
Fall 1 =2 beschrankt. 

In §9, B und C behandeln wir zwei Iterationsverfahren fiir endliche Ring- 
gebiete, wovon das eine (KomAaTU) schon in [5] auf Konvergenz untersucht war. 
Wegen der iiberaus raschen Konvergenz sei auf das Verfahren in § 9, C besonders 
hingewiesen. 

$6. Spiegelung von Gebreten 


Gegeben seien zwei beliebige fremde Gebiete G, und G, der z-Ebene, die eine 
gemeinsame Randkomponente C haben sollen. Ist C ein Kreis, so ist klar, was 
darunter zu verstehen ist, wenn wir sagen, G, und G, seien beziiglich C spiegel- 
bildlich. Wir erklaren die spiegelbildliche Lage von G, und G, beziiglich C all- 
gemein wie folgt. 
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Definition 1. Die Gebiete G, und G, mit der gemeinsamen Randkomponente C 
heiBen dann spiegelbildlich beztighich C, wenn es eine konforme Abbildung w =f (2) 
von GyrCuG, gibt, so daB 

(a) f(C) ein Kreis K der w-Ebene ist, und 

(B) #(G,) und f(G,) beztighch IK spiegelbildlich lregen. 

Die beziiglich K spiegelbildliche Lage der Bilder wird also durch f+ auf die 
spiegelbildliche Lage der Urbilder iibertragen. 

Manchmal ist es notwendig zu wissen, ob G, iiberhaupt ein Spiegelbild hat, 
d.h. ob G, an C spiegelungsfihig ist. Aus unserer Definition geht unmittelbar 
hervor, daB dies genau dann der Fall ist, wenn es eine konforme Abbildung 
w=f(z) von G, auf ein Gebiet G, gibt derart, daB 

(x) #(C) ein Kreis K ist, und 

(6) die Umkehrfunktion f4(w) in das beziiglich A genommene Spiegelbild G3 
von Gi hinein analytisch und schlicht fortsetzbar ist. Dann ist /*(G:) =G, das 
beziiglich C genommene Spiegelbild von G,. 


Wir bemerken sogleich mehrere einfache Folgerungen aus der Definition. 


(a) Damit G, an C spiegelungsfahig ist, ist notwendig, daB C eine analytische 
Jordan-Kurve ist. 


(b) Die durch f1Sf/ angegebene ,,natiirliche’ Zuordnung spiegelbildlicher 
Punkte ist indirekt konform. Dabei bezeichnete S die Spiegelung am Kreis K, 
die ja indirekt konform ist. 


(c) Insbesondere ist G3 indirekt konformes Bild von G,. Ist daher speziell G, 
ein zweifach zusammenhangendes Gebiet vom Modul M(>1), so ist auch G, 
zweifach zusammenhangend und vom selben Modul. 


(d) Liegen G, und G, beziiglich C spiegelbildlich, und wird das Gesamtgebiet 
G,VCuUG, einer konformen Abbildung / unterworfen, so liegen h(G,) und h(G,) 
beziiglich h(C) spiegelbildlich. Denn fh bildet (G,) und h(G,) auf Gebiete ab, 
die beziiglich K spiegelbildlich liegen. 


(e) Beztiglich einer Randkomponente C von G, existiert héchstens ein Spiegel- 
bild. Bezeichnen wir namlich (ftir den Fall zweier Spiegelbilder) mit F(z) und 
Fy(z) (2€G,) die in (b) genannten Zuordnungen spiegelbildlicher Punkte, so ist 
F, (2) =4(2) fiir z€G, aus einer hinreichend kleinen Umgebung von C, da wir 
ja eine Spiegelung an einer analytischen Kurve vor uns haben, welche eindeutig 
erklart ist ([7], S. 360). Also ist die in G, analytische Funktion [A @) —Be)| =o 
fiir alle hinreichend nahe an C liegenden z€G,. Daher ist F (z) a ig (2 € G,) 
und folglich sind die Spiegelbilder beziiglich C eindeutig festgelegt. y, 


(f) Wir betrachten ein an C spiegelungsfahiges Gebiet G,; sein Spiegelbild 
sel Gy. Wird nun G, allein einer konformen Abbildung w =h(z) unterworfen 
so ist h(z) genau dann tiber C hinaus nach G, schlicht fortsetzbar, wenn h(Gy) 
an h(C) spiegelungsfahig ist. : 

(1) Ist A(z) nach G, schlicht fortsetzbar, so wende man (d) an. 
(ii) Sei h(G,) an h(C) spiegelungsfahig, und * bezeichne die Spiegelung an C 
bzw. an A(C). Dann betrachten wir fiir z€G, die Funktion 
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die in G, analytisch ist. Bei Annaherung an C ist g(z) stetig (C und /(C) sind 
analytische Jordan-Kurven), und es ist g(z) =h(z) (¢€C). Nach dem Stetig- 
keitsprinzip setzt g(z) die gegebene Funktion h(z) nach G, hinein fort, und es 
ist h(G,) gleich dem Spiegelbild von h(G,) beziiglich h(C). 

Spater benétigen wir vor allem folgende Teilaussage von (f): 

(g) Die Funktion w=h(z) bilde G, konform auf Gj ab. G, sei an C und Gj 
an h(C) spiegelungsfahig. Dann ]a8t sich h(z) in das Spiegelbild G, von G, be- 
ziiglich C hinein analytisch fortsetzen, und es ist 4(G,) das Spiegelbild von h (G,) 
beziiglich h(C). 

Bei unseren Untersuchungen in § 7 bis § 9 wird der Fall von Bedeutung sein, 
daB ein Gebiet G mehrfach spiegelungsfahig ist. 


Definition 2. Das n-fach zusammenhangende Gebiet G heiPt einmal allseitig 
spregelungsfalig, wenn es beztiglich jeder seiner n Randkomponenten spiegelungs- 
fahig ist. Ist dies der Fall, und ist auBerdem jedes der n gewonnenen Spiegelbilder 
von G selbst wieder einmal allseitig spiegelungsfahig, so nennen wir G zweimal 
allseitig spiegelungsfahig. 

Hohere Spiegelungsfahigkeit wird entsprechend erklart. Durch mehrfache 
Anwendung von (g) erhalt man allgemeiner 

(g’) Die Funktion w =h(z) bilde G, konform auf G, ab. G, und Gy seien jeweils 
k-mal allsettig spiegelungsfalig; die durch thre Spiegelbilder erweiterten Gebiete 
heiBen G, und G3. Dann laBt sich h(z) nach G, hinein analytisch fortsetzen und 
es ist h(G,) =Gs: 


§7. Das ,,tterterende Verfahren“ von KOEBE 


A. Angabe des Verfahrens. Vorgegeben sei ein beliebiges m-fach zusammen- 
hangendes Gebiet G=Gy, der %-Ebene mit co€G») und den Randkomponenten 
OM (¢=1, 2,...,), welches auf ein Vollkreisgebiet abzubilden sei. Bei dem 
ersten hier zu behandelnden Verfahren geht man dazu folgendermaBen vor. 


Beim ersten Schritt werde das gesamte AuBere von C(° (also ein einfach 


zusammenhangendes Gebiet) durch z,=/, (2%) konform und normiert (um Pyeaco 
Pilt. hy(Z) =z +2 2 a) auf das AuBere eines Kreises CY abgebildet; 
die Bilder pe-ciennen oe mit h,(G)) =G, und h,(C®) =C. Das Gebiet G, ist 
bei nicht-spezieller Wahl der C® (was auch im folgenden angenommen wird) 
genau einmal spiegelungsfahig, namlich am Kreis C{”. 

Beim zweiten Schritt werde das gesamte AuBere von C$) durch z,=h,(%) 
konform und normiert auf das AuBere eines Kreises CY’ abgebildet; wir setzen 
hy(G,) =Gy und hy(C) =C®. Die Kurve C{? verliert bei diesem Schritt zwar 
ihre Kreisgestalt, jedoch bleibt (§6, (d)) G, auch tiber CP hinaus spiegelungs- 
fahig. Neu hinzu kommt die Spiegelungsfahigkeit von G, tiber CY). Und so fort. 

Beim (mn +1)-ten Schritt betrachten wir wieder die ,,erste“ Kurve C”) und 
bilden ihr gesamtes AuBeres durch z, ,,;=,+1(%,) konform und normiert auf das 
AuBere eines Kreises C{"t” ab; wir setzen h,,,,(G,) =G,+, und h,.,(C”) =Cet. 

Wie das Verfahren fortfahrt, ist klar; die entstehenden Bilder der m Ausgangs- 
kurven werden periodisch durchlaufen und stets das gesamte AuBere eines solchen 
Bildes konform und normiert auf das AuBere eines Kreises abgebildet. Nach 


dA 
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allgemein m Schritten sind m konforme Abbildungen /, (29), 22 (%), 73 (22); voy 
Ay (Zn—1) Gurchgefiihrt, und durch Hintereinanderschalten derselben erhalten wir 


die m-te Iteration: 
fn (Zo) = Pom «++ Re My (Zo) (29 € Go; Wire ply oa ah 


Diese Funktion bildet G) konform und normiert auf G,, ab, und wir werden in 
Teil B sehen, daB f,,(z)) unsere gesuchte Abbildungsfunktion von Gp auf das 
Vollkreisgebiet mit wachsendem m immer besser approximiert. ; 
Zuvor wollen wir jedoch untersuchen, wie oft G,, allseitig spiegelungsfahig 
ist. Dazu betrachten wir die 7-te Kurve, die G,,_, berandet, also Core (A475). 
G,,_, sei tiber sie, d.h. wenn man nur Spiegelbilder von G,,_, innerhalb Cre 
betrachtet, S’”~-mal allseitig spiegelungsfahig. Bei unserem m-ten Schritt 
werde nun C!"—» auf einen Kreis abgebildet. Diese Abbildung laBt die Spiege- 
lungen von She an C{"-) (i+-4,) bestehen (§6, (d)), d.h. G,, ist SY) = S%"-»- 
mal iiber C\”) hinein allseitig spiegelungsfahig (7+:7)). Uber den Kreis By ist 
G,, dagegen ersichtlich S$?” =4+ min S(-)_mal allseitig spiegelungsfahig. Be- 


zeichnet also insgesamt S’”) die allseitige Spiegelungsfahigkeit von G,, tiber C”, 
und wird beim m-ten Schritt Cae in einen Kreis tibergefiihrt, so besteht die 
Rekursion 

Ses ae AS150; 144) 


v 


(7-4) 


Unser Gebiet G,, ist dann genau (min S%”)-mal allseitig spiegelungsfahig. Aus 
(7.4) findet man nun leicht, daB G,, fiir m =n erstmals einmal allseitig spiegelungs- 
fahig ist. Fiir m=2mn—1 erhalt man erstmals zweimalige allseitige Spiegelungs- 
fahigkeit, und allgemein ist G,, fiir m—=k-n—(k—1) erstmals k-mal allseitig 

spiegelungsfahig. Anders ausgedriickt: 
Satz 7.1. Das nach m Iterationsschritten aus Gy hervorgegangene Gebiet Gy, 
ist N(m) = = 
n— 


: [nal allseitig spiegelungsfahig. 


B. Konvergenz des Verfahrens. Nachdem wir nunmehr untersucht haben, wie 
die allseitige Spiegelungsfahigkeit von G,, mit m zunimmt, sind wir in der Lage, 
daraus Schliisse zu ziehen auf die Annaherung unserer gesuchten Normalabbil- 
dung durch },, (2). 

Hierzu beweisen wir folgenden allgemeinen Satz ®, 


Satz 7.2. Das Vollkreisgebiet V der w-Ebene sei durch z=g(w) konform und 
normiert auf das n-fach zusammenhdngende Gebiet G(co €G) der z-Ebene abgebildet. 
Die Rénder C, bzw. C; von V bzw. G liegen in |w|SR bzw. |z|SR, und G set 
N-mal allseitig spiegelungsfihig. Weiter sei D =dist(V; UC) gesetzt, wenn CM 


die durch enmalige Spiegelung von V an seinen Randkreisen entstehenden neuen 
Kreise sind. Die Kreise C,; seien ferner durch n Jordan-Kurven I; umschlossen, 


so dap C; und I; n fremde Ringgebiete von den Moduln M;(>1) beranden, und 
schlieBlich ser M = min M,. 


1Si<n 


6 Die in Satz 7.2 und dessen Beweis verwendeten Bezeichnungen decken sich nicht 
mit denen, die in der Hauptentwicklung von § 7 vorkommen. 
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Dann gilt mit der Konstanten K(M) aus Satz 1.2 


. 2 2 
7.2) |g(w) ~w| =o HARROD IPL wer). 
Der Satz prazisiert, wie sich hohe Spiegelungsfahigkeit von G niederschlagt 
in einer guten Approximation der Identitat w durch die Abbildungsfunktion g (w). 
Man sieht auBerdem, daB die Auswirkung hoher Spiegelungsfahigkeit von G auf 
gute Approximation von w durch g(w) um so gréRer ist, je groBer M ist, das 
heiBt (in einem relativen Sinne) je weiter die Kurven C,; voneinander weg liegen. 


Fig. 2 


Beweis. Wir fithren zunachst einige Bezeichnungen ein. Nach Vorausset- 
zung, bzw. weil V Vollkreisgebiet ist, sind G und V N-mal allseitig spiegelungs- 
fahig. Die durch ihre Spiegelbilder erweiterten Gebiete nennen wir G* und V*, 
und § 6 (g’) ist zu entnehmen, dafB-g(w) in V* hinein schlicht fortsetzbar ist und 
als Bild g(V*) =G* liefert. Die nach einmaliger allseitiger Spiegelung von V 
bzw. G an den Kurven C;=C bzw. C;=C® entstehenden n(n —1) neuen Kur- 
ven nennen wir C® bzw. CY. Eine weitere allseitige Spiegelung an jenen Kurven 
ergibt die Rander C®) und C. SchlieBlich wird G* von den Kurven C®) und V* 
von den Kreisen C\%) berandet. 

Weiter setzen wir g(I}) =, und verfolgen nun die I] und J} beim Spiegelungs- 
prozeB. Spiegelung von J} = 2 an C;=C liefert y™ ; alle weiteren ersten Spiegel- 
bilder der I} (wenn also J} an C; gespiegelt wird (i=+j)) nennen wir I}. Die 
n(n —1) Kurven I umschlieBen je ein Cf und beranden mit ihm ein Ring- 
gebiet vom Modul = M, wogegen ein y" stets m—1 Kurven J)” umschlieBt und 
mit C!° ein Ringgebiet vom Modul = M berandet. Der weitere SpiegelungsprozeB 
fiihrt entsprechend zu den J{®) und y’ (vy =1, 2,...,.N), und in der z-Ebene er- 
halten wir ganz analog NM? =¢(I}) und y®=g¢(y®) (v=1, 2, ..., N). 

Nun gilt nach der Cauchyschen Integralformel fiir jedes we V 


esl g(o)—@ 1 g(a)—-@ 
heh 271 i; o—w AE 271 fs Aide 


uc) |o|=e 
a a 
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wobei @ hinreichend gro8 ist. Wegen der Normierung von g(w) in w=oo ba 
schwindet das zweite Integral, und da w auBerhalb der Integrationskreise CH 
liegt, kénnen wir weiter schreiben 


pigs cathe &(o)—Ai gg 
g(w) —w = oe pa f ag ae 
v cm 


mit era Konstanten A,, die wir jetzt geeignet wahlen. Wenn wir namlich 
mit d\) den Durchmesser von C$ C™) bezeichnen, so gibt es nach dem schon in § 5 


dN) 
verwendeten Satz von Junc ([13], S.140) einen Kreis mit Radius —— ick vero 


enthalt; A, sei sein Mittelpunkt. Wir setzen noch /)= Lange von C{®) und er- 
halten damit fir wEV 


(73) |g) —w| Se 2M) GM) ss ae Dia) sta Oa 


weil ja |w—@|=D gilt (wEV, wEC§). 
a) Beurteilung von > (%))2. Bezeichnen wir im folgenden mit F(C) stets die 


V 
von C umschlossene Flache, so kénnen wir unter Verwendung des Satzes von 
CARLEMAN (Satz 1.1) sagen: 


Dr An CMs 4 


vee = M2 M4 
v 4 
1 
<e 
SS SAN = aI, aie F(T) 
1 
a ya ee Fy (0) 42 ; 
40 pe cara anda Sa DEG )s Man TR; 


man beachte, daB y®) und [}’) (y=N,N—1,...,2) stets ein Ringgebiet vom 
Modul = M2 beranden. 
f) Beurteilung von Dd )?. Hier miissen wir auBer dem Satz von CARLEMAN 


noch den Satz 1.2 herancielions sonst ist die Abschatzung ganz analog zu «). 
Wir erhalten 


AN) 2< a 2, ah? 
2a") = K(M) DF(R™) S---S K(M)-M. 


Insgesamt ergibt sich mit «) und ) aus (7.3) 


: {40-4+K(M) M®} R? 
|g(w) —w|<+* aE j aN (wEV), 


womit Satz 7.2 bewiesen ist?. 


7 Auf demselben Weg kann man auch zum Beispiel | g’(w) — 1| abschatzen. 


Durchfiihrung der konformen Abbildung 107 


Nunmehr sind alle Vorbereitungen getroffen, um die Konvergenzgeschwindig- 
keit des Koebeschen iterierenden Verfahrens zu beurteilen®. 


Satz 7.3. Gegeben sei ein beliebiges n-fach zusammenhingendes Gebiet G der 
2-Ebene (co€ G) mit den Randern C,, die alle in \z|<R liegen mégen. Jedes Rand- 
kontinuum C; sei durch eine Jordan-Kurve I; umschlossen, so dab C,; und Ty 


n fremde Roepe von den Moduln M;,(>1) beranden; wir setzen M = min M,. 


1Sisn 
Wetter bezeichne w=f(z) die konforme und in z=co normierte Abbildung von 


G auf ein Vollkreisgebiet V, und es sei D =dist (V, U y;) gesetzt, wenn y, die durch 
eimmalige Spiegelung von V an dessen Randkreisen entstehenden n(n —1) neuen 
Kreise sind. 

Dann gilt fiir die durch das Koebesche iterierende Verfahren gewonnenen Néahe- 
rungsfunktionen f,,(z) die Abschatzung 


Kil 2. R2 
eee (a) = F(ayy ae a ge a ey et Dg a 
3 ut laa 
K(M) bedeutet wieder die Konstante aus Satz 1.2. 


Beweis. Zu beurteilen ist 


[fm (2) — f(2)| =|fnf*(@)— | =|g(@)—w| ~~ we), 


m—| (Satz 7.1). Als die Kurven J; 


7 


Minds Satz 7.2. 1st anwendbar mit N=| 


in Satz 7.2 hat man dabei jetzt f(I}) zu nehmen. Weiter ist noch zu beachten, 
daB alle Rander von V und von G,,=},,(G) in |w|<2R bzw. |z,,| 2K liegen, 
was wegen der Normierung von f(z) und f,,(z) in zoo der Fall ist ([10], S. 178). 


$8. Das ,,Iterationsverfahren“ von KOEBE 


Neben dem in § 7 besprochenen Verfahren hat KOEBE 1920 noch ein weiteres 
funktionentheoretisches Iterationsverfahren zur konformen Abbildung mehrfach 
zusammenhiangender Gebiete eingefiihrt, dem wir uns jetzt zuwenden wollen. 


A. Angabe des Verfahrens. Vorgegeben sei ein beliebiges m-fach zusammen- 
hangendes Gebiet G=G, der z-Ebene mit co€ Gy und den Randkomponenten 
C® (¢=1, 2,...,), welches auf ein Vollkreisgebiet abzubilden sei. 

Beim ersten Schritt werde das gesamte AuBere von C{) durch z,=h, (z) kon- 
form und normiert auf das AuBere eines Kreises K abgebildet und sodann das 
entstehende Gebiet G,=/,(G)) an K gespiegelt. Der Kreis K% wird geléscht, 
so daB als Ea a se iibrig bleiben Ci) =h, (Cora) Re erases nie und 
deren Spiegelbilder: C2, ,,; als Spiegelbild von C!) an K® (i =1,2,...,m—1). 

Allgemein liegen vor dem m-ten Schritt die Kurven C{”"~ GA, 2. Area) 
vor, und durch z,,=h,,(%,—1) Wird nun das gesamte er von C{” Nyaa ete 
und normiert auf es eee eines Kreises K pa und sodann das ent- 
stehende, von K) und h,,(C%"~)) (i =2,3, ..-,%m—1) berandete Gebiet an K 
gespiegelt. Der Kreis K‘ (nm) yy geloscht, so jap Ak Berandungen tibrig Deen 
CM—h, (CM?) (¢=1, 2, ...,%,-1—1) sowie deren Spiegelbilder an To eeine. 
gesamt also 7,,=27,-1— 2 Randor 


8 Wir verlassen jetzt die Bezeichnungen von Satz 7.2 wieder. 
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Wie in §7 erhalt man durch Hintereinanderschalten der m Abbildungen 
hy (Zo), Hg (21) +++» Ym (Zn—1) die m-te Naherung fiir die gesuchte Abbildungsfunk- 
tion 

fon (Zo) = Pm + + + Mig My (2p) (29:6 Gog eps BG 


wir setzen noch G,,=f,, (Go). 


B. Konvergenz des Verfahrens. Wie in § 7 ist auch hier fur die Konvergenz- 
geschwindigkeit die allseitige Spiegelungsfahigkeit N(m) von G,, maBgebend, die 
wir daher zunachst untersuchen. 

Dazu legen wir jeder Berandung Cc’ eine Kennzahl bei, die angibt, wie oft 
G,, gespiegelt werden mu8, um C () 7a erreichen. 

Nach m—1 Schritten liegen vor C{"-?,...,C@"7. Wir nehmen an, daB 
ihre Kennzahlen mit wachsendem unterem Index monoton (im weiteren Sinne) 
zanehmen; dies ist bei nicht-spezieller Wahl von G, richtig fiir 7 =1. Von diesen 
Randern haben 

A, die Kennzahl k 
A,,, die Kennzahl k + 1 


A die Kennzahl s. 


s 


Beim m-ten Schritt wird zunachst C’"~) auf den Kreis K) abgebildet, und 
es verbleiben 
A,—1 Rander mit der Kennzahl k 


(8.1) A,,, Rander mit der Kennzahl k + 1 
A. Rander mit der Kennzahl s. 


Die Cie) ah, (rey) (¢ =1, 2,...,7%,-1—1) haben namlich dieselbe Kennzahl wie 
(Sete (§ 6, (d)), also wachsen auch deren Kennzahlen mit unterem Index mono- 
ton. Durch die Spiegelung dieser Kurven am Kreis K) entstehen 


A,—1 Rander mit der Kennzahl (2k +1) + 
(8.2) A,,, Rander mit der Kennzahl (2k+1)+k+41 

As Rander mit der Kennzahl (2k +1)+s. 
a die Kennzahlen der Spiegelbilder wachsen monoton, so daB insgesamt die 
C™ ((=1, 2,...,7%,) mit unterem Index monoton wachsende Kennzahlen be- 
sitzen. 

Bei diesem Proze8 wird also auf systematische Art stets eine Randkomponente 
kleinster Kennzahl auf einen Kreis abgebildet und durch dessen Léschung weg- 
geschafft. Durch (8.1) und (8.2) ist ein rekursives Verfahren angegeben, mit 
dessen Hilfe sich berechnen la8t, wieviel Rinder nach dem m-ten Schritt eine 
gewisse Kennzahl haben. Das Minimum aller auftretenden Kennzahlen ist die 


allseitige Spiegelungsfahigkeit N(m) von G,,, und da beim Beweis von Satz 7 


nur die Spiegelungsfahigkeit von G,,, nicht aber die spezielle Konstruktion des 
Iterationsverfahrens bentitzt wurde, kénnen wir sagen: 
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Satz 8.1. Unter den Bedingungen von Satz 7.3 gilt die dortige Aussage (7.4), 


—1 es A 4 . 
wenn ee | durch die zum Koebeschen ,,Iterationsverfahren‘ gehorige Funktion 


N(m) ersetzt wird. 


Nun ist in dem hier vorliegenden Fall kein handlicher Ausdruck fiir N(m) 
angebbar, auBer im Fall m =2, auf den wir in § 9 zuriickkommen werden. Fiir 
n> 2 beschranken wir uns auf die Feststellung, daB N(m) im allgemeinen be- 
a4 in §7. Ist G) zum Beispiel dreifach 
zusammenhangend, so ergeben namlich unsere Rekursionsformeln (8.1) und (8.2), 
daB 


trachtlich langsamer wachst als 


N(m) = 12 2B 4 See OMe Loy Oe e110 
in §8 erstmals fiir m= Op Se ee Ne - Bile Sy EA Ph OE 
in §7 erstmals ftir m= Be 5 7 Oni H1oeMsat 19 21 


gilt. Dieser Nachteil des Iterationsverfahrens riihrt daher, daB8 durch die Spie- 
gelung an K™) zwar die Kennzahlen der Rander vergréBert werden, jedoch auch 
die Anzahl der wegzuschaffenden Rander exponentiell wachst, wahrend beim 
iterierenden Verfahren stets nur m Rander in Betracht bleiben. Dadurch wird 
dort die allseitige Spiegelungsfahigkeit gleichmaBiger erhoht als hier. 

Wir konnen also feststellen, daB vom praktischen Standpunkt aus das 
,iterierende Verfahren“ (§ 7) dem ,,Iterationsverfahren“ im Fall n=3 vorzu- 
ziehen ist, wenngleich auch jetzt N(m) 7 co (m-—>oo) gilt und damit die Kon- 
vergenz des Verfahrens feststeht. 


$9. Iterationsverfahren fiir Ringgebiete 


Wir kommen jetzt speziell zur konformen Abbildung zweifach zusammen- 
hangender Gebiete (= Ringgebiete). In A sollen die Ergebnisse von § 7 und § 8 
auf den Fall » =2 spezialisiert werden, wahrend wir in B und C die konforme 
Abbildung von Ringgebieten auf konzentrische Kreisringe behandeln. 


A. Der Sonderfall n=2 in § 7 und § 8. Vorgegeben sei ein beliebiges zweifach 
zusammenhangendes Gebiet G der z-Ebene (co€ G) mit den Randern C, und C,, 
die beide in |z|<R liegen mogen. Der Modul von G sei M(>1), und das Voll- 
kreisgebiet, auf das G vermége w =/(z) konform und in z = oo normiert abgebildet 
werden soll, heiBe V. Mit D bezeichnen wir wieder den Abstand von V zu den 
beiden Kreisen, die durch einmalige Spiegelung von V an seinen beiden Rand- 
kreisen entstehen. 

Um f(z) zu approximieren, kénnen wir zwei Methoden verwenden. 

a) Zunachst wird das iterierende Verfahren von §7 jetzt zu einem alter- 
nierenden Verfahren, und die danach gebildeten Naherungsfunktionen fiir /(z) 
bezeichnen wir mit f,,(z). Um Satz 7.3 anwenden zu kénnen, bemerken wir, 
daB sich zwischen C, und C, eine Niveaulinie I" legen 1aBt, derart, daB C, und I’ 
sowie C, und J" zwei fremde Ringgebiete des Moduls /M beranden; |/M tritt 
daher an die Stelle von M in § 7. Wir erhalten somit 
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Satz 9.1. Bildet man fiir das Gebiet G des Moduls M die Funktionen f,,(z) nach 
dem iterierenden Verfahren von § 7, so gilt 


0.4) [fat -1@)| SURO my | (EG: mB). 


Die Abschatzung lieBe sich noch etwas verbessern, da wir nur die allseitige 
Spiegelungsfahigkeit ausgeniitzt haben und daher eine weitere Spiegelungsmég- 
lichkeit unberiicksichtigt lieBen. 

b) Nun wenden wir auf G das Iterationsverfahren des §8 an. Wie schon 
dort bemerkt, werden jetzt durch die bei jedem Schritt eingeschaltete Spiegelung 
an einem Kreis die Kennzahlen der nicht geléschten Rander erhéht. Da nun 
im Fall 2 =2 dieser Gewinn nicht mit einer Vermehrung der zu loschenden Rander 
verbunden ist, ergibt sich eine betrachtliche Konvergenzverbesserung. 

Wir ermitteln zunachst N(m) mit Hilfe unserer Rekursionsformeln (8.1) und 
(8.2). Vor dem m-ten Schritt habe der eine Rand die Kennzahl «,,_,, der andere 
die Kennzahl «,,_»<«,,_,- Dann haben die Rander nach dem m-ten Schritt die 
Kennzahlen «1 bZW. o,=2 %m—2+%,—-1+1. Nennen wir also die Kennzahl des 
beim m-ten Schritt neu hinzukommenden Randes  (m), so gilt p (m) =2q(m—2)+ 
y(m—1)+1; als Anfangsbedingungen dieser Rekursion haben wir (1) =1, 
y(2)=2. Die Lésung der Differenzengleichung ergibt sich zu 


OU oe ae as 


und da N(m) =g(m—1) ist, haben wir damit 


m+1 Sa \ ee 
(9.2) NG ae 
gefunden. Wir bemerken noch 
2h Phi 
N(m) 2—— ~ 26 (m = 2) 


Ganz entsprechend wie in a) erhalten wir somit 


Satz 9.2. Bildet man fiir das Gebiet G des Moduls M die Funktionen Pata) 
nach dem Iterationsverfahren von § 8, so gilt 


(9.3)  |fn(2) —#(2)|S Gee pep EGS MED), 


wober N(m) aus (9.2) zu entnehmen ist und augerdem N(m) 22 (m= 2) gilt. 


Eine geringe Mehrarbeit beim Iterationsverfahren (zu der bei jedem Schritt 


auszuftihrenden konformen Abbildung kommt noch eine Spiegelung) bewirkt also 
eine erhebliche Beschleunigung der Konvergenz. 


B. Das Iterationsverfahren von KOMATU. Fiir den Rest dieses Paragraphen 
ist unser Ausgangsgebiet G=G, der z-Ebene ein Ringgebiet, das auBen von 
|2o0| =1 und innen von einem in |zo|<1 verlaufenden Kontinuum Cy berandet 
ist; Z%=O sel Gy. Ist M(>1) der Modul von Go, so liegt die Aufgabe vor, Go 


durch eine konforme Abbildung w=f/(z 9) auf den konzentrischen Kreisring 
M+1<|w|<1 so abzubilden, daB f(4) =4 wird. 
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Als eine erste Moglichkeit zur Lésung dieser Aufgabe besprechen wir das 
Verfahren von Komatu (siehe [5]), das auf derselben Grundidee beruht wie das 
iterierende Verfahren fiir 7 =2. Der erste Iterationsschritt wird in zwei Etappen 
ausgefiihrt. Zunachst sei 4=g,(%) eine konforme Abbildung des gesamten 
AuBeren von Cy auf |¢,|>4 mit g,(co) =co. Sodann bilde man durch 2 =h, (4) 
das gesamte Innere von g,(|z)| =1) auf |z,|<4 ab mit /,(0) =0; auBerdem sei 
h, so normiert, daB h,g,(1) =1 wird. Die analytische Jordan-Kurve h,({f,| =1) 
nennen wir C,; mit |z,|=1 zusammen berandet sie G,=/,g,(G,). 

Allgemein wird beim m-ten Schritt zunachst das ganze AuBere von C,,_, 
durch t,,=8m(Zm—1) auf |¢,,|>1 abgebildet (g,,(00) =0c), und hernach tibe man 
eine konforme Abbildung z,,=,,(f,,) aus, die das Innere von g,,(|z,,-1|=1) auf 
l2m|<4 abbildet, wobei h,,(0) =0 und A,,g,,(1) =1 sein soll. Als Bild von |¢,,| =1 
erscheint C,,, die mit |z,,| =1 zusammen G,, berandet. 

Durch Hintereinanderschalten der 2m ersten Abbildungen erhalten wir die 
m-te Iteration 


Tin (20) = Pm Gin +» Ne So My &y (2) (2 € Gy; m21). 


In einer fritheren Arbeit ([5]) hatten wir durch wiederholte Anwendung einer 
Verfeinerung des Schwarzschen Lemmas gezeigt, daB die Kurven C,, in einem 
prazisen Sinne immer mehr kreisformig werden, woraus sich dann die Konvergenz- 
aussage 


(9.4) \/m (20) — f@)| =O({* arctg Mall") — (z4€ Gy; moo) 


ergab. Die folgenden Uberlegungen verbessern diese Abschatzung. 


Man stellt namlich durch eine Betrachtung der beiden Einzelschritte fest: 
Wenn G,, iiber C,,, hinein «,,-mal spiegelungsfahig ist, so ist G,,,, tiber C,,,, hinein 
(,,+-2)-mal spiegelungsfahig. Das bedeutet, daB G,, iiber C,, hinein 2m-mal 
gespiegelt werden kann, oder, daB die durch z,,=/,,/+(w) vermittelte Abbildung 
von M1<|w|<1 auf G,, schlicht nach (M7)?”<|w|<1 fortgesetzt werden kann 
(§6, (g)). Dabei erkennt man weiter, daB bei dieser Fortsetzung z,,=0 ausge- 
lassen wird. Wir kénnen daher Satz 5.1 auf die Funktion /,, 4 (w) und fir7 =M~?” 
anwenden und erhalten so 


inf 2 (w) ae ww | = C sae? (Mt< | w| <1; m= 1,2, as 
das heiBt 
NCA are heey orate (n6Gg 7 = 2). 


Satz 9.3. Es sei G ein Ringgebiet der z-Ebene, das augen von |z|=1 und 
innen von einem in |z|<1 verlaufenden Kontinuum berandet sev; z=O sev ¢G. 
Bildet man dann nach dem Verfahren von Komatu die Iterationen f,,(2) fiir die 
konforme und durch f (1) =1 normierte Abbildung w =f (z) von G auf M+<|w|<1, 
so gilt 
(9.5) | fm (2) — f(z)|SC-M-*™ @EG3 M==172,49.),, 


wobei C<12,6 die in Satz 5.1 auftretende absolute Konstante ist?. 


® Beim Vergleich mit Satz 9.1 beachte man jedoch, daf jetzt ftir einen Iterations- 
schritt zwei konforme Abbildungen auszufiihren sind. 
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Unsere neue Konvergenzaussage (9.5) ist besser als die frihere (9.4), da 


Mt< + arc tg M+ ist. Wir bemerken noch, daB ftir die Giltigkeit von (9.5) 


PLUS Tis (1) —{ wichtig war, nicht jedoch die Erfiilltheit aller vorhergegangenen 
Normierungen /, g;, (41) =1. 

C. Ein weiteres Iterationsverfahren. Wir bemerkten bereits, daB das Koebe- 
sche Iterationsverfahren dem iterierenden Verfahren im Falle der Ringgebiete tiber- 
legen ist. Da nun das Verfahren von Komatu in volliger Analogie zum iterierenden 
Verfahren aufgebaut ist, entsteht die Frage, ob man nicht in naheliegender Weise 
durch Einschalten einer Spiegelung bei jedem Einzelschritt des Komatuschen 
Verfahrens ebenso eine Konvergenzverbesserung erzielen kann?°. Dies ist der 
Fall, wie wir jetzt naher ausfiihren wollen. 

Unser Ausgangsgebiet sei wie in B beschaffen, doch spiegeln wir gleich Co 
an |z)|=1 und erhalten Cy. Jeder Iterationsschritt zerfallt nun wieder in zwei 
Teilschritte. Zunadchst bilde man das ganze AuBere von Cy durch 4,=g, (29) 
(g,(co) =o) konform auf |f,|>4 ab und spiegele g,(Co) an |4|=1, was Co er- 
gebe. Danach werde das ganze Innere von g,(Co) durch z,=/,(4) konform auf 
|z,|<1 abgebildet (A, (0) =0), wobei Co in C, tibergehe; das Spiegelbild von C, 
beziiglich |z,|=1 heiBe Cj. Beim nachsten Schritt tibernehmen C, und Cj; die 
Rolle von C, und Cg, und das Verfahren fahrt in Analogie zum Iterationsver- 
fahren fort. Nach m Iterationsschritten (jeder bestehend aus zwei Teilschritten) 
erhalten wir durch Zusammensetzen unserer konformen Abbildungen 


Fi, (Zo) = Him 8m + x 82h & (Zo) (Z€Go; m =4,2,...), 


und durch F,, (29)/F, (4) ist auch noch der Normierungsvorschrift (z9=1) + (z,,= 1) 
genugt. 

Wir vermuten, da8 in Analogie zu A.b schon diese Funktionen F,, (z9)/F,, (1) 
mit wachsendem m die gesuchte Abbildung w=f(z)) von Gy auf M1<|w|<1 
immer besser approximieren, doch haben wir dafiir nicht die nétigen Hilfsmittel 
bereit. Daher schalten wir nach Ausfiihrung der durch F,,(z)) ausgedriickten m 
Iterationsschritte noch eine weitere konforme Abbildung @,,(z,,) vor, die das 
ganze Innere des AuBenrandes von F,,(G)) auf den Einheitskreis abbildet, und 
unsere m-te Iteration heiBt nunmehr endgiiltig 


fn (Zo) = Pry Min Sm ++ Ire SoMa 81 (20) (%€ Go; m = 4, 2,...); 


die Normierung von @,, soll so gewahlt werden, daB PD, (0} =0' und: f2(1) =4 
wird. Zur Gewinnung von f,,(z 9) sind also 2m-—+1 konforme Abbildungen und 
2m Spiegelungen notwendig. 

Die Einschaltung dieser Spiegelungen hat nun zur Folge, daB unser gewon- 
nenes Bildgebiet G,,=},,(Gp) tiber seinen Innenrand hinein 


(9.6) Nim) =2(4"—1)  (m=1) 


mal gespiegelt werden kann. Diese Formel fiir N(m) findet man aus einer Re- 
kursionsformel [N(m + 2) =5 N(m +1) —4N(m)], oder auf Grund der Bemerkung, 
daB sich hier nach m Schritten dasselbe N ergeben muB wie in A (9.2) fiir 2m +1. 


Ee ; 
Die Untersuchung eines solchen Verfahrens wurde durch Herrn HUBNER an- 


, gerege. 
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Wir erhalten daher genau so wie in Teil B den 


Satz 9.4. Es sei G ein Ringgebiet der z-Ebene, das aupen von |z| =1 und innen 
von einem in |z|<1 verlaufenden Kontinuum berandet sei; z=0 sei €G. Bildet 
man dann nach dem oben geschilderten Verfahren die Iterationen f,,(z) fiir die kon- 
forme und durch {(1)=1 normierte Abbildung w= f(z) von G auf M41<|w|<1, 
so gilt 


(9.7) lfn(2) —f(2)| S$C-MN™ CCG mm AD Sy, 


wobet C<12,6 die in Satz 5.1 auftretende absolute Konstante und N(m) aus (9.6) 
zu entnehmen ist. 


Die Einschaltung von Spiegelungen bei jedem Iterationsschritt wurde also 
mit einer betrachtlichen Beschleunigung der Konvergenz des Verfahrens belohnt. 


§ 10. Anhang 


Im folgenden sollen die in § 9, B und C geschilderten Iterationsverfahren noch 
auf das Beispiel eines exzentrischen Kreisrings angewandt werden. In diesem 
Falle lassen sich die exakte Abbildung und die Naherungsfunktionen und damit 
auch der Approximationsfehler explizit berechnen. Wir vergleichen die Ergeb- 
nisse der verschiedenen Iterationsverfahren und machen insbesondere Aussagen 
iiber die Qualitat der Fehlerabschatzungen (9.5) und (9.7) [vgl. die Zusdtze zu 
den Satzen 9.3 und 9.4, unten]. 


A. Das Verfahren von KOMATU. a) Aufstellung der Naherungen. Unser Aus- 
gangsgebiet G=G, der z-Ebene sei als Bild von M4<|w|<1 vermége 


1—@y, (Wao 
1—a@) 1—G)w 


(10.1) = = |, (w) 

dargestellt. Es wird von |z)| =1 und einem z)=0 umschlieBenden Kreis berandet, 
falls noch | a)|<M~ gewdahlt ist. Die nach dem Komatuschen Verfahren (§ 9, B) 
konstruierten Gebiete G,, sind wieder exzentrische Kreisringe, deren eine Beran- 
dung |z,,| =1 ist, so daB G,, als Bild von M4<|w|<1 vermége 


_1-—a@, WG, 
1—Qy, 1—Ay,W 


(10.2) Bs = 1, (w) 


angesetzt werden kann, und dabei ist gleich w=1 in z,,=—1 tibergefithrt. Zur 
Beurteilung des Verfahrens ermitteln wir {a,,} rekursiv. 


Zunachst wird die Abbildung z,,—>t,,,, in den beiden Schritten z,,—>¢,,,1 und 
bn41—tm41 Aurchgefiihrt, wobei 


, af 7 ie sh eeele 
bim+1 = M1, (Zp) , Lat = 


gilt. Um z,,= 0c in t,,,,;= 0c zu bringen, ist offenbar 


gi 1 wat! Gm 
fi43(00) Migi(ooc) MM’ 


am 


alSOmmog— ar 


erforderlich. Die Abbildung Zy)44=/yy41(tm+41) ist wieder linear, mit /,,4(0) =0. 
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Nun 1aBt sich a,,,, leicht errechnen: 


, Ay , a a. 
em+1 — 0 ag ee =0 > bad ry ar > on = by (tin41/4) 77 Ln (An) Se ae aa 
: F a 
Andererseits ist auch w=a,,,, Urbild von z,,,;=0, d.h. es gilt Omti=Fqe oder 
allgemein 
(10.3) a,= 2 (wm = 0;1,2,...)— “mit |a)|)=— uM 


Unsere Naherungsfunktionen f,,(%) fiir die exakte Abbildungsfunktion 
} (2) =l9* (%) schreiben sich dann in der Form f,, (Z9) —1,J¢*(z)) und sind durch 
(10.1) bis (10.3) vollstandig bestimmt; f,,(1)=1 ist ersichtlich erfillt. 


b) Beurteilung der Approximation an einem konkreten Beisprel. Zu unter- 
suchen ist der Maximalfehler 


by, = a3 | f(a) — fo] = max [2 — | = MAX ln () — 3 
dabei kénnen wir, da /,,(w) —w sogar in |w|<1 regular ist und daher nach dem 


Maximumprinzip max | #,,(w) —w| mit ymonoton wachst, auch nur a L,,(w) —w| 
w|=r w|= 


ermitteln. Ist a) und daher a,, reell (zur reellen Achse symmetrischer Kreisring), 
so ergibt eine leichte Rechnung 


2a 
Mem 


0, 24, = 


Ist unser Ausgangsgebiet Gy) etwa das von |z)|=1 und |z)+34| = 3{ beran- 
dete Gebiet (M = 2, ay= 4), so ergibt sich als maximale Abweichung der Komatu- 
schen Naherung gegen die exakte Abbildungsfunktion 


(10.4) Oa Ao (= O° ARO 
c) Beurteilung dey Giite der Fehlerabschitzung (9.5). Wir untersuchen nun 
= a |/,,(w) —w| im allgemeinen Fall eines komplexen a. Es gilt 


1—GQy, W—Gy, 
1—G,, 1—G6,,wW 


=|(1—4,,) (1+ 4,,) (w — a,,) (1 +4,,w) — w| + O(M-4™) 
= | (Zp @ + ay) (w — 1)| + O(M-*™) 
= Fo) et Nira 

bm = 4|ao| -M-?™ 4+ O(M-4™) (it Seen als 


[An (2) — w| = 


WwW 


dabei wurde (10.3) verwendet und 4)=1|a9| gesetzt. Da |ay| nur der Bedingung 
|@)|<M™ unterworfen ist und M>1 beliebig sein kann, sehen wir an diesem 
Beispiel: 

Zusatz zu Satz 9.3. Die Grigenordnung der A bschatzung (9.5) ist bestméglich. 


Weiter kann die dort auftretende absolute Konstante C = 12,6 durch keine Konstante 
<4 ersetzt werden. 
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B. Das Verfahren von §9, C. a) Aufstellung der Naherungen. Als Ausgangs- 
gebiet G=G, wahlen wir wieder das Bild von M1<|w|<1 vermége 


Ei1=a ee Oy ne Me ; 4 
(10.5) ay (ag ate 1p (w) ; | 4o| <M. 
Unser erstes Ziel ist es, die Funktion F,,(z9) zu bestimmen (Bezeichnungen wie in 
§9,C). Diese Funktion bildet ganz allgemein Gp auf ein Ringgebiet ab, das 
ber seinen Innenrand N(m) mal (bis zur Kurve C,,) und iiber seinen AuBenrand, 
bis zu |z,,|=4 hin betrachtet, K(m) mal spiegelungsfahig ist, wobei 


(10.6) K(m) =3(4"—1) und N(m) =2(4"%—1) (i SOFAS 2a.) 


gilt. Fiihrt man die entsprechenden Spiegelungen in der w-Ebene durch, so ent- 
stehen die beiden Kreise |w| =7,, und |w| =R,,, als Urbilder von C,, und |z,,| =1, 
und dabei ist 


(10.7) (pe Sind ie, Via 


Setzen wir daher ¢,,=R,, z,,, So ist der Zusammenhang zwischen ¢,, und w von 


der Form 
Cea Patinlty = 1,,(w). 


” Ry =a w 


Hierin ist R,, durch (10.7) bekannt, so daB noch a,, zu ermitteln ist. 
Dazu fithren wir den (m-+1)-ten Iterationsschritt aus, dessen ersten Teil- 
schritt z,,—>t,,,, Wir in den beiden Etappen z,,>¢,,41 und ¢),41>ty)4, erledigen. 


Wir setzen ise er 
, ey he =a m 5 
Lm +1 =i lin (Gy) ) Ln+1 a IC ) 


C,, wird dadurch in |¢,,,,;|=1 tbergefiihrt und z,,= oo in 4,,,,;= 00, wenn 


= 1 Ym Ym o> Vm 
= = ee) Oe 5 ACL LO, = a, 
tm +1 (0) Ea((ES)) eee eines 


gewahlt wird. Die Abbildung z,,.,;=/,,13 (4,41) ist wieder linear, mit h,,,,(0) =0. 
Damit laBt sich a,,,, errechnen: 


, Vin x. 2 
em+1 = k)' > Ent Sah inti Go A R. “Ay, =O, ae Lin (Qin tae) p 
m 


Andererseits ist auch w=a,,.,R,,1, und also €,,=1,(@y41R,y41) Urbild von 
240, —0, also gilt 


2 

Vm 

a ate Ay OGY Gygaq -_@,,, 
aT Ugg 2 ea Rien 


und (10.6) und (10.7) liefern sofort 
See ie eh OUC Bi ig Vigna) CoE eee 
Unsere gesuchte Funktion F(z) stellt sich somit dar durch 
FE, (20) = Rin mn bo * (20) » 


und die auftretenden GréBen R,, und a,, sind durch a) und M bestimmt. Aus 
dieser Funktion stellen wir nun, gemaB unserem Vorgehen in § 9, C, zweierlei 
Naherungen her. 
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Erstens wollen wir lediglich noch der Normierungsvorschrift (%)=1) > (m=) 

geniigen, indem wir 

Q(z.) — F,(zq)/F,(1) = =). (m = 0,1,2,--.) 

In (Zo) We MNS OME 1 c Ip 3(1) 
betrachten. Die zweite Naherung entsteht durch Vorschalten der Funktion DP, 
mit ®,,(0) =0, die das ganze Innere des AuBenrandes von £;, (Go) auf den Ein- 
heitskreis abbildet, wobei ®,,F,,(1) =1 werden soll. Dies leistet aber 
at ee 
1—A  41—As’ 


Deaton |= Boonie (8) = A=L,*(0) = (Ol eudons 


und unsere zweite Naherungsfunktion lautet somit in unserem Fall 
ie (2) a ®,, By (Zo) = Lin Ig * (2) % 


b) Beurteilung der Approximation an einem konkreten Beisprel. Wie in A 
betrachten wir den Sonderfall, daB Gy von |z)| =1 und |z)+ 24 | = 34 berandet ist 
(M =2, ay=}). Als Maximalfehler 


Oy max Lea) — f ()| = max |L,, lo * (Zo) — 40 * (Zo) ee) = 


ZyE Go 2 EG 
erhalten wir sofort 
(10.9) Of == 2 Ace 20 Ree ted ia (a 0, Aer, ad 


Fiir den Maximalfehler von /® (z,) gegen /(z)) bestimmen wir 


(1) =i) = Lyn (w) | 
max | ft,’ (2) — lo (2 m w| = max 
ZE Go i/ ( o y ( 0 wl =1 | Zn (1) |o|=1 


= |a)| ¢,,- max |w — 1| + O(é,) 
|w|=1 


Rn—G4n em om Rin 
1—Ay Ry Ry—GAynw 


= 2 hag| Em + O(e.,), 


mit ¢,,=M—~3@"-)_ Dabei braucht a, nicht reell zu sein. Hieraus folgt nun 
fiir unser Beispiel (M =2, ay=4) 


(10.10) Of =o dE Ua APO See 


c) Beurteilung der Giite der Fehlerabschdtzung (9.7). Lassen wir zur Definition 
von Gp rein imaginar Werte von dy zu, so finden wir wie in A allgemein 


Om = 2|ao|- M804 O Ae) 


, 


52 = 4] a] -M-8¢"-0 4 O(M-#@"-2), 


Insbesondere ergibt die Abschatzung von 6) wie oben den folgenden 


Zusatzzu Satz 9.4. Die Gripenordnung der A bschatzung (9.7) ist bestméglich. 
Weiter kann die dort auftretende absolute Konstante CX12,6 durch keine Kon- 
stante <4 ersetzt werden. 


Unsere Untersuchung hat auch ergeben, da8 im Falle des exzentrischen Kreis- 


rings das Vorschalten von @,, iiberfliissig war; die ohne diese Funktion gewonnenen 


Naherungen /\}) (z)) lieferten sogar eine bessere Approximation von f(z) als {® CAE 
falls a) rein imaginar war. 
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Ein Vergleich der geschilderten Iterationsverfahren und der Ergebnisse (10.4), 
(10.9), (10.10) liefert insgesamt folgendes Bild: 


—_eeeee—ee————— ss eee 


Maximaler Approximationsfehler 
im Sonderfall M=2; a)= } 


Anzahl der bis zum m-ten Schritt 
notwendigen Operationen 


Komatu-Verfahren 


2m konforme Abbildungen 


Verfahren von § 9, C 


ohne Verwendung von ®,, |~™” 


| 2m konforme Abbildungen 


2m Spiegelungen 


Verfahren von § 9, C 


2m -+ 1 konforme Abbildungen 


: 62d gabe) 
mit Verwendung von @,, Es 2m Spiegelungen 
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